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Abstract 


Approximate  Theories  of  Elastic  Rods  With  Apphcations 

by 

Jeffrey  Stephen  Turcotte 
Doctor  of  Philosophy  in  Mechanical  Engineering 
University  of  California,  Berkeley 
Assistant  Professor  Oliver  M.  O’Reilly,  Chair 

In  this  dissertation,  approximate  theories  involving  combinations  of  small  and 
moderate  strains  and  rotations  for  elastic  rods  are  developed.  Their  usefulness  is 
illustrated  with  several  applications.  The  rod  theory  used  to  construct  these  theories 
is  the  directed  (or  Cosserat)  rod  theory  developed  by  Green,  Naghdi  and  several  of 
their  co-workers.  The  approximate  theories  which  aie  developed  in  this  dissertation 
are  rendered  properly  invariant  under  arbitrary  superposed  rigid  body  motions  by 
extending  some  recent  work  of  Casey  and  Naghdi,  and  of  O’Reilly.  These  extensions 
were  developed  to  render  the  properly  invariant  theories  more  amenable  to  applica¬ 
tions. 

The  approximate  theory  that  is  the  primary  focus  of  this  work  is  one  involving 
small  strain  and  moderate  rotation.  A  parallel  development  for  a  directed  surface 
was  performed  earher  by  Naghdi  and  Vongsarnpigoon.  Specifically,  it  is  shown  that 
there  are  considerable  simphfications  in  the  balance  and  constitutive  laws  as  well  as 
in  the  strain-displacement  relations  because  of  the  assumptions  made  in  this  theory. 

Among  the  apphcations  considered  is  a  discussion  of  certain  flexural  modes  of 
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vibration  which  have  escaped  attention  in  the  literature  and  which  are  also  present 
in  a  Timoshenko  beam.  Of  additional  interest  is  the  small  strain-moderate  rotation 
theory,  and  three  specific  examples  of  its  use  are  presented.  These  applications  show 
the  nature  of  the  nonlinearity  introduced  by  moderate  rotation  as  weU  as  some  restric¬ 
tions  of  the  theory  which  are  dependent  on  the  specific  application.  This  dissertation 
concludes  with  a  discussion  of  the  whirling  rod.  Specifically,  this  problem  is  formu¬ 
lated  as  one  of  infinitesimal  deformations  superposed  on  a  large,  steady  deformation. 
Among  the  issues  addressed  is  the  influence  of  the  steady  deformation  on  the  small 
amplitude  vibrational  response  of  the  rod. 


Assistant  Professor  Oliver  O’Reilly,  Chairman 
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Chapter  1 


Introduction 

1.1  General  Background 

A  rod  theory  is  a  theory  in  which  the  field  equations  have  been  reduced  to  dependence 
on  a  single  spatial  coordinate  {e.g.,  the  arc  length).  Primarily,  there  have  been 
two  approaches  to  the  development  of  rod  theories  in  the  literature.  One  approach 
starts  with  three-dimensional  continuum  mechanics  and  is  based  on  an  approximation 
procedure.  Examples  of  various  procedures  may  be  found  in  Love  [41,  §251  -  §254] 
and  in  Green,  Naghdi  and  Wenner  [25].  The  second  approach  is  referred  to  as  the 
direct  approach,  where  the  balance  laws  are  postulated  a  priori^  and  the  concept  of  a 
directed  curve  is  used.  According  to  Ericksen  and  TruesdeU  [15],  directed  media  were 
first  considered  by  Duhem  [13]  and  later  studied  by  E.  and  F.  Cosserat  [8].  The  theory 
of  a  directed  or  Cosserat  curve  has  been  extended  by  A.  E.  Green,  P.  M.  Naghdi  and 
others.  A  review  of  this  development  is  provided  by  Naghdi  [46].  The  theory  used  in 
this  dissertation  has  two  directors  attached  to  a  space  curve,  although  theories  with 
an  arbitrary  number  of  directors  are  developed  in  [20]  and  [25]. 
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1.2  Scope  and  Content  of  this  Dissertation 


In  this  dissertation,  we  are  primarily  concerned  with  developing  approximate  theories 
of  Cosserat  curves,  but  we  also  provide  results  for  a  number  of  applications.  Approx¬ 
imate  rod  theories  that  are  based  on  three-dimensional  continuum  mechanics  have 
been  addressed  by  a  number  of  authors  such  as  Shield  and  Im  [57]  and  Danielson 
and  Hodges  [10].  These  theories  typically  allow  for  large  or  moderate  rotations  but 
require  the  strain  to  remain  small,  thereby  encompassing  more  deformations  than  the 
lineax  theory  at  the  expense  of  adding  some  complexity  or  nonlinearity.  The  approx¬ 
imate  theories  of  a  Cosserat  curve  we  develop  in  this  dissertation,  however,  are  quite 
distinct  from  other  theories.  We  give  special  attention  to  the  approximate  theory  of 
small  strain  accompanied  by  moderate  rotation,  as  the  infinitesimal  theory  has  been 
previously  detailed  in  [18],  [20],  [22]  and  [26],  and  theories  involving  moderate  strain 
lack  sufficient  development  of  their  constitutive  responses  to  be  useful  at  this  time. 

In  Chapter  2,  we  recall  all  of  the  field  equations  of  the  purely  mechanical  two- 
director  model  of  a  constrained  Cosserat  curve.  We  also  discuss  the  correspon¬ 
dence  between  the  direct  approach  and  three-dimensional  continuum  mechanics.  In 
addition,  we  discuss  invariance  requirements  and  the  properly  invariant  theory  of 
O’Reilly  [50]  for  approximate  theories  of  an  elastic  rod.  In  the  process,  we  modify 
an  assumption  in  [50]  that  is  overly  restrictive  on  the  invariance  requirements  for  the 
constraint  forces. 

Chapter  3  details  four  approximate  theories  of  Cosserat  curves  after  deriving  the 
differential  equation  for  the  rotation  tensor  of  the  Cosserat  curve  in  terms  of  the  strain 
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measures.  The  four  theories  result  from  all  combinations  of  small  and  moderate 
strain  and  rotation.  In  this  chapter,  which  is  motivated  by  the  work  of  Naghdi 
and  Vongsaxnpigoon  [49],  the  differential  equation  for  the  rotation  tensor  is  used  to 
form  the  approximation  for  the  cases  having  small  strain.  This  method  proves  to  be 
unsatisfactory  in  the  approximations  involving  moderate  strain,  and  leads  us  to  use 
an  alternate  approach  which  we  discuss  in  Chapter  4. 

Following  the  work  of  Casey  and  Naghdi  [5]  in  three-dimensional  continuum  me¬ 
chanics,  we  use  the  symmetric  and  skew  symmetric  parts  of  a  displacement  gradient  to 
represent  strain  and  rotation  in  Chapter  4.  This  provides  a  straight-forwaxd  method 
of  developing  the  same  four  approximate  theories  we  identified  in  Chapter  3.  We 
then  show  that  this  method  yields  the  same  results  as  the  method  of  Chapter  3  to 
the  order  of  approximation.  We  also  use  the  method  to  show  that  these  approximate 
theories  are  not  properly  invariant  unless  we  use  the  auxiliary  motion  of  O’Reilly  [50], 
which  we  summarize  in  Chapter  2. 

In  Chapter  5,  we  develop  a  constrained  theory  to  further  simplify  the  balance 
laws.  We  then  proceed  to  specialize  to  the  case  of  straight  rods.  This  allows  for 
significant  simplification  of  the  balance  and  constitutive  laws.  For  future  reference, 
the  field  equations  of  the  infinitesimal  theory  are  also  recorded  in  this  chapter. 

Chapter  6  deals  with  hnear  extensional  and  flexural  vibrations.  Since  the  linear 
flexural  balance  laws  of  a  Cosserat  curve  with  two  directors  are  equivalent  to  the 
Timoshenko  beam  theory,  much  of  the  development  given  for  free  flexural  vibrations 
is  available  in  the  hterature.  However,  we  do  determine  and  discuss  certain  modes  that 
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have  not  been  previously  explored,  and  the  example  given  serves  as  a  preparation  for 
another  example  in  Chapter  7.  The  free  extensional  response  that  we  obtain  solutions 
for  is  unique  in  its  inclusion  of  non-axisymmetric  lateral  extensions.  In  fact,  we  use  a 
rectangular  cross-section  in  our  example  and  can  see  clear  phase  relationships  between 
the  longitudinal  and  lateral  extensions. 

We  provide  examples  of  the  moderate  rotation  theory  in  Chapter  7.  Two  of  these 
examples  highlight  the  coupling  between  flexural  and  extensional  response,  which  is 
one  of  the  primary  differences  between  an  infinitesimal  theory  and  a  moderate  rotation 
theory.  We  provide  both  a  static  example  and  a  free  vibration  example  and  complete 
the  chapter  with  an  example  showing  the  sensitivity  of  the  theory  to  deformations. 
This  last  example  shows  that  the  theory  is  not  valid  for  all  deformations  of  a  specific 
rod. 

In  Chapter  8  we  develop  the  balance  laws  for  a  rod,  one  end  of  which  rotates  at  a 
constant  angular  velocity,  and  consider  the  free  vibrations  superposed  on  the  steady 
deformation.  After  a  brief  review  of  the  literature  on  this  subject,  we  formulate  the 
appropriate  boundary  and  initial  value  problem  for  this  system.  We  then  establish  the 
balance  laws  for  the  associated  steady,  large  motion.  Finally,  we  obtain  the  balance 
laws  for  the  superposed  vibrations.  The  resulting  equations  are  rendered  properly 
invariant  using  modified  auxiliary  motions.  We  do  not  solve  for  the  vibration  modes 
of  the  rod  because  the  exphcit  form  of  the  free  energy  is  not  assumed.  However,  we 
discuss  an  approximate  solution  procedure  that  could  be  used  if  this  function  were 
available. 
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1.3  Notation 


In  this  dissertation  we  use  several  standard  notational  conventions,  some  of  which 
are  summarized  here.  All  lower  case  Latin  indices  range  from  1  to  3,  lower  case 
Greek  indices  range  from  1  to  2,  upper  case  Latin  indices  range  from  1  to  R,  and  the 
summation  convention  for  repeated  indices  is  employed.  The  tensor  product  a  0  b  of 
any  two  vectors  a  and  b  is  a  second  order  tensor  defined  by  the  operation  (a  0  b)c  = 
(c  •  b)a  for  every  vector  c.  The  Euclidean  norm  of  a  second  order  tensor,  denoted 
II  •  II,  is  defined  by  ||  A  |p=  <r(AA^)  =  A*  A,  where  “tr”  is  the  trace  operation,  and 
indicates  the  inner  product  of  two  tensors.  If  e  =  e{t)  =  sup^  II  ^  ll?  where  “sup” 
stands  for  the  supremum  (or  least  upper  bound)  of  a  non-empty  bounded  set  of  real 
numbers,  then,  following  Casey  and  Naghdi  [6],  we  use  the  notation  h(x)  =  0(e") 
as  e  — >  0  for  a  function  h(x)  which  is  defined  in  the  neighborhood  of  x  =  0  if  there 
exists  a  constant  C  >  0  such  that  ||  h(x)  ||<  Ce".  For  further  details  on  notation  and 
inequalities,  we  refer  the  reader  to  [6,  §2.1]  and  Naghdi  and  Vongsarnpigoon  [48,  §4]. 
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Chapter  2 


Theory  of  a  Cosserat  or  Directed  Curve 

In  this  chapter  we  summarize  the  theory  of  a  Cosserat  curve  developed  by  A.  E. 
Green,  P.  M.  Naghdi  and  their  co-workers.  This  theory  was  developed  in  several 
papers;  however  the  main  results  can  be  found  in  [17],  [19],  [20],  [26]  and  [46]. 

2.1  Kinematics 

We  recall,  from  Naghdi  [46],  the  concept  of  a  Cosserat  curve.  This  curve  consists  of  a 
material  curve,  which  is  embedded  in  Euclidean  three-space  and  to  which  at  each 
material  point  of  the  curve  a  set  of  deformable  vector  fields  or  directors  is  defined. 
For  purposes  of  the  present  discussion,  it  suffices  to  consider  the  case  where  there  are 
just  two  directors  di  and  d2.  We  employ  a  convected  coordinate  i  to  identify  the 
points  of  the  material  curve. 

The  motion  of  the  Cosserat  curve  is  defined  by  the  vector- valued  functions 

r  =  r(^,t),  d„  =  d,(^,i),  (2.1) 

where  a  =  1  and  2,  r(^,t)  uniquely  identifies  the  present  configuration  I  of  the 
material  curve,  and  (2.1)  uniquely  identifies  the  present  configuration  of  the  Cosserat 
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curve.  It  is  eissumed  that  the  scalar  triple  product 


[dij  d2)  da]  >  Oj 


(2.2) 


where  a  standard  notational  convention  is  adopted: 


(2.3) 


We  recall  that  the  assunaption  (2.2)  implies  that  a  set  of  reciprocal  vectors  {d‘}  can 
be  uniquely  defined: 

d*'  •  d,-  =  8),  (2.4) 

where  8*^  is  the  Kronecker  delta.  The  material  time  derivatives  of  the  quantities  r  and 
da  are 

V  =  f(^,t),  Wa=da(6<),  (2.5) 


where  the  superposed  dot  denotes  the  partial  derivatives  of  these  functions  with  re¬ 
spect  to  t  (keeping  ^  fixed). 

For  convenience,  a  fixed  reference  configuration  of  the  Cosserat  curve  is  also  de¬ 
fined  by  the  vector- valued  functions 


R  =  R(0,  Da  =  Da(0- 


(2.6) 


It  is  assumed  that 


where,  paralleling  (2.3), 


[Di, D2, D3]  >  0, 


D3  = 


98,  ‘ 


(2.7) 


(2.8) 
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Again,  (2.7)  implies  that  a  set  of  reciprocal  vectors  {D'}  can  be  defined; 


D‘  •  Dj-  =  S'j.  (2.9) 

We  now  recall  the  kinematic  tensors  of  the  Cosserat  curve  which  were  introduced 
by  Naghdi  [46,  §13]:^ 


F  =  d„  ®  D“  +  da  0  D^, 

(2.10) 

C  =  F’’F,  E=i(C-I), 

(2,11) 

oGa=  _  0D^ 

(2.12) 

K«  =  F^G„  -  oG„. 

(2.13) 

In  (2.10)  and  (2.12),  (g)  is  the  standard  tensor  product  of  two  vectors.  In  addition,  we 
recall  the  definition  of  the  displacement  vector  of  the  material  curve  i: 

n{C,t)  =  T{U)-m)-  (2-14) 

With  the  assistance  of  (2.3)  and  (2.8)  it  follows  that 

^  =  da  -  Da.  (2.15) 

By  substituting  (2.10)  into  (2.11)  -  (2.13),  we  also  obtain 

E=-  (di  •  dj  -  D,-  •  Bj)  D‘  D^'  =  EijD*  ®  (2.16) 

K“  =  ^  ‘  ®  ® 

^Naghdi  [46]  does  not  explicitly  introduce  the  tensor  E  in  direct  notation  and  the  symbols  ■yij/2 
for  the  components  of  this  tensor  are  usually  used  in  the  literature.  The  notation  E  was  introduced 
by  O’Reilly  [50]. 
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2.2  Balance  Laws 


Prior  to  discussing  the  balance  laws  of  a  Cosserat  curve,  the  folloAving  five  fields  are 
introduced:  n,  the  contact  force;  k",  the  intrinsic  director  forces;  and  m",  the  director 
forces.  We  also  define  A  =  A(0  =  PoiO'/^  ^  length  of  the  rod 

in  the  reference  configuration,  where  D33  =  D3  •  D3. 

From  Nagh-di  [46,  §9]  and  Green  and  Naghdi  [20],  the  balance  laws  for  a  Cosserat 
curve  are,  mass  conservation: 

A  =  0,  y“  =  0,  =  (2.18) 


balance  of  linear  momentum: 


—  +  Af  =  A  (v  +  y  w„) , 


(2.19) 


(two)  balances  of  director  momentum: 


5m“ 


+  Al"  -  k"  =  A  (y“v  +  , 


(2.20) 


and  the  balance  of  moment  of  momentum: 


da  X  n  +  d„  X  k“  +  ^  X  m“  =  0.  (2.21) 

In  (2.18)  -  (2.21),  f  =  f(^,t)  is  the  assigned  force,  1“  =  l“(^,t)  are  the  assigned 
director  forces,  and  y"  =  y“(0  and  y“^  =  y^°‘  =  y“^(0  are  inertia  coefficients. 

To  obtain  one  component  form  of  the  balance  laws,  following  Green,  Naghdi  and 
Wenner  [26],  the  various  vectors  are  resolved  onto  {d^}:  e.g.,  n  =  n'd,-,  k"  =  ^"*d, 
and  in“  =  m“*d,.  For  convenience,  we  define  the  differences  between  the  assigned 
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forces  and  acceleration  terms  as 


f  =  f  -  (v  +  y"w„)  =  /’di, 

(2.22) 

q"  =  1“  -  (i/"v  +  =  q“‘d,-. 

(2.23) 

The  restdting  component  forms  of  (2.19)  -  (2.21)  are 

Sn*  _ 

^  +  A^n*-  +  A/‘  =  0, 

(2.24) 

r.  + 

(2.25) 

^  -  m“^A„^  =  0, 

(2.26) 

jbV  _  pA  ^  =  0, 

(2.27) 

where  A/  =  d’  •  ddrjd^. 


2.3  Constitutive  Equations  and  Constraints 

In  this  section  we  outhne  the  constitutive  and  constraint  responses  of  a  Cosserat 
curve.  This  theory  is  given  in  Naghdi  [46,  §10]  and  was  further  developed  by  Naghdi 
and  Rubin  [47]  and  O’Reilly  [50].  We  begin  by  postulating  the  existence  of  a  free 
energy  il)  of  the  Cosserat  curve  and  R  mechanical  constraints  having  the  following 
functional  dependences;^ 

1*  =  ^  =^(F,G<,.„G.,D„e).  (2.28) 

=v''(F,G„„G.,D,,{)  =  0.  (2.29) 

^See  [47]  for  examples  of  such  constraints. 


10 


As  V’  is  a  constitutive  response  for  the  Cosserat  curve,  it  is  assumed  to  be  objective: 

V)'*'  =  ■0(QF,  QG„,  oGa,Dj,^)  =  i^  (F,  Gc,  oG^jD,-,^) ,  (2.30) 

where  Q  =  Q{t)  is  a  proper  orthogonal  tensor  representing  the  rotation  due  to  the 
superposed  rigid  body  motion,  and  the  transformation  of  F  and  G^  under  superposed 
rigid  body  motions  is  explained  fully  by  O’Reilly  [50,  §2].  O’Reilly  [50]  and  others  have 
also  shown  that  the  objectivity  of  tp  can  be  used  to  reduce  the  functional  dependence 
of  V’  to: 

^  {Eij,  Kaj,  oGa,  Di,  0  »  (2-31) 

where  Eij  and  Kai  =  «a.3  are  the  respective  components  of  E  and  as  defined  in 
(2.16)  and  (2.17). 

Recalling  the  expression  for  mechanical  power  (c/.  Naghdi  [45,  §10]): 

n  •  ^  +  k"  •  w„  +  m“  •  =  A^,  (2.32) 

dC 

and,  if  we  assume  that  the  kinetical  responses  may  be  separated  into  their  constitutive 
responses  and  constraint  responses,  respectively,  as 

n  =  n  +  n,  k“  =  k“  +  k“,  m“  =  m"  +  m",  (2.33) 


then  the  constitutive  responses,  obtained  using  a  procedure  similaj  to  that  discussed 
by  Naghdi  and  Rubin  [47],  are  given  by 


n  =  A 


(  d\l)  dxj) 

\dE^k  "  dU' 


(2.34) 


k"  =  A 


/  .  dtp  dd/jA 

^  dK^0  )  ’ 


(2.35) 
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(2.36) 


dib 

ih“  =  A-^ - djfc. 

OKak 


Following  [50]  and  [24],  we  assume  that  each  of  the  constraint  functions  is  objective. 
Using  standard  arguments,  they  must  also  have  the  representation 


oGaj  DtiO 


(2.37) 


and  the  constraint  responses  are 


ra  (  I  -  „  aoL 

A  r. 

m  =Pl^ - clfc=Pi,C  ) 

OKak 


(2.38) 

(2.39) 

(2.40) 


where  pL  =  PL(^,t)  Lagrange  multipliers  and  we  have  used  the  summation  con¬ 
vention  for  the  index  L  {L  —  1, . . . ,  R). 

The  component  forms  of  the  constitutive  equations  for  isothermal  rods  are  given 
in  Green  and  Naghdi  [20,  §13]  and  in  Green,  Laws  and  Naghdi  [19,  §8]  as 

dij) 


A  3  A  a3  \  3  \ 

n  -m  A„,  - 


2  dEaz  ’ 

=  2X 

dijp 


m°'*  =  A 


dKai  ’ 


(2.41) 

(2.42) 

(2.43) 

(2.44) 


where  the  function  ip  must  be  written  in  a  form  which  allows  for  the  appropriate 
symmetries  of  Eap.^  Clearly,  the  component  form  of  the  constraint  responses  can  be 
presented  in  a  similar  manner. 


®We  show  later  an  acceptable  form  for  ip. 
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2.4  Correspondence  with  Three-Dimensional 


Continuum  Mechanics 

In  this  section,  we  show  how  the  Cosserat  curve  may  be  used  to  model  a  three- 
dimensional  continuum  B.  Let  be  a  set  of  convected  coordinates  for  B.  These 
coordinates  uniquely  define  the  material  points  of  B.  The  position  vector  p  =  p{d*,t) 
of  a  material  point  of  B  is  assumed  to  be  approximated  by 

p  =  r  -f  O^da,  (2-45) 

where  0^  =  ^.  A  similar  situation  holds  for  the  reference  configuration  (but  there  the 
corresponding  approximation  is  exact;  i.e.,  it  is  a  representation). 

To  help  the  reader  gain  insight  into  the  theory  and  to  provide  formulae  which 
are  useful  in  applications,  we  recall  here  all  of  the  basic  kinetical  integrations  from 
the  three-dimensional  theory.  For  this  purpose,  we  assume  as  in  Green,  Naghdi  and 
Wenner  [25]  that  the  region  occupied  by  the  body  at  time  t  is  a  neighborhood  of  £ 
bounded  by  a  material  surface  F{6^,  ^^,0  —  0.  Then  ^  =  constant  defines  a  material 
plane  A  which  is  bounded  by  F  =  0.  Following  [25,  §2],  the  curve  £  is  fixed  relative 
to  the  surface  F  =  0  by  the  conditions 

where  g*  is  the  determinant  of  the  metric  tensor  g*j  =  gj  •  gj  for  the  three-dimensional 
convected  coordinate  basis  {g»}  =  P*  is  the  mass  density  of  the  body,  and 

y°‘  =  Since  p*^/^  is  independent  of  time  due  to  the  conservation  of  mass. 
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the  curve  £  is  a  material  curve.  Using  these  definitions,  the  kinetical  vectors  of  the 
Cosserat  curve  theory  can  be  calculated.  The  resulting  expressions  are  again  taken 
from  [25].  The  contact  force  is  given  by 


n  = 


(2.47) 


the  intrinsic  director  forces  are  given  by 


k“  =  J 


(2.48) 


and  the  director  forces  are  given  by 


m 


a 


d9^de\ 


(2.49) 


where  T*  =  the  stress  vectors  associated  with  the  contravariant  com¬ 

ponents  T*-’  of  the  Cauchy  stress  tensor  T.  The  assigned  forces  axe,  again  according 
to  [25,  §2], 


Af  =  I l^p*y/^rd9U9^  +  f  [(t'  -  A'T")  d9^  -  (t^  -  A^T")  d9^] ,  (2.50) 

Xr  =  I  j^p*y/rr9'^d9U9^+ j  9°‘  [(t^  -  A^T^)  d9'^  -  (t^  -  A^'T^)  d9^]  ,  (2.51) 


where  f*  is  the  three-dimensional  body  force,  and  the  line  integrals  are  taken  along 
the  curve  defined  by  the  intersection  of  the  surfaces  (  =  constant  and  F  =  0.  Also, 
A  =  A“ga  -t-  g3  is  a  vector  which  is  tangent  to  the  surface  F  =  0.  Finally,  the  inertia 
coefficients  are  specified  (according  to  Green  and  Naghdi  [20,  §10])  by  (2.46)  and 

A  =  y  p*y/^  9^9^  d9U9^.  (2.52) 
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2.5  Invariant  Theory 


O’Reilly  [50]  has  shown  that  the  infinitesimal  theory  of  an  elastic  Cosserat  curve  is  not 
properly  invariant  under  arbitrary  superposed  rigid  body  motions.  The  approximate 
theories  to  be  considered  in  this  work,  which  include,  but  are  not  limited  to,  the 
infinitesimal  theory,  suffer  from  the  same  drawback  (we  will  show  this  to  be  the  case 
in  section  4.5).  Following  the  work  of  Casey  and  Naghdi  [5],  which  was  in  the  context 
of  three-dimensional  continuum  mechtinics,  O’Reilly  [50]  defines  an  auxdhary  motion 
of  a  Cosserat  curve  which  is  properly  invariant  and  yet  able  to  incorporate  large 
motions.  Since  we  will  use  this  theory,  as  well  as  a  slightly  modified  form  of  it,  we 
provide  an  outhne  of  the  theory  here  and  refer  the  reader  to  [50]  for  additional  details. 
We  also  discuss  a  modification  to  the  theory  given  in  [50]  in  which  the  invariance  of 
the  Lagrange  multipliers  in  a  constrained  theory  is  not  imposed.  These  invariance 
requirements  are  discussed  in  detail  by  O’ReiUy  and  Turcotte  [52]. 

2.5.1  Invariance  Requirements 

We  begin  with  a  discussion  of  the  invariance  requirements.  For  an  unconstrained  rod 
under  superposed  rigid  body  motions,  following  the  assumptions  of  Naghdi  [45,  §8] 
for  the  case  of  a  Cosserat  surface  and  Green  and  Laws  [17]  for  the  case  of  a  Cosserat 
curve,  O’ReiUy  [50]  assumes  the  kinetical  fields  transform  under  superposed  rigid 
body  motions  according  to 

n+  =  Qn,  k“+  =  Qk“,  =  Qm“.  (2.53) 
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Thus,  these  fields  are  termed  objective.  However,  in  the  case  of  a  constrained  rod,  as 
shown  by  O’Reilly  and  Turcotte  [52],^  it  is  not  appropriate  to  require  the  full  kinetical 
fields  n,  k",  m“  to  be  objective.  The  parts  of  these  fields  determined  by  constitutive 
equations  are  seen  to  be  objective  by  first  recognizing  that 

df  =  Qd,,  (2.54) 

which  implies  that  the  measures  E  and  are  invariant.  Since  we  have  assumed  that 
the  free  energy  V’  is  also  invariant,  the  partial  derivatives  dtl^ldEij  and  drpIdKai  are 
completely  invariant.®  Using  this  result  and  (2.54)  in  the  constitutive  equations  (2.34) 
-  (2.36)  shows  that  the  constitutive  equations  are  objective: 

n+  =  Qn,  k“+  =  Qk“,  m“+  =  Qih“,  (2.55) 

but  the  Lagrange  multipHers  which  determine  the  constraint  response  need  not  be 
objective  (remain  constant  under  superposed  rigid  motions),  depending  upon  how  the 
superposed  motion  is  brought  about. 

For  example,  if  the  superposed  rigid  body  motions  are  caused  by  the  addition 
of  (obviously  non-objective)  surface  tractions  which,  by  careful  construction,  do  not 
affect  the  constitutive  responses,  the  Lagrange  multipHers  will  generally  not  be  ob¬ 
jective.  To  show  this  clearly,  consider  the  balance  laws  for  two  motions  that  differ  by 
a  superposed  rigid  body  motion: 

^This  point  was  first  made  by  Casey  and  Carroll  [4]  in  the  context  of  three-dimensional  continuum 
mechanics. 

®For  completeness,  we  should  also  note  that  the  measures  A^!  are  unaltered  under  superposed 
rigid  body  motions. 
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(2.57) 


^  ^  +  Al»  -  it”  -  k“  =  A  (!,“v  +  y‘l>wg) 

^  ^  +  Af+  =  A  (v+  +  y-wi)  ,  (2.58) 

fl^  +  +  Al"+  -  k“+  -  t"*  =  A  (y<-v*  +  .  (2.59) 

Application  of  the  invaiiance  requirements  (2.55)  and  substituting  for  dn/d^  from 
(2.56)  into  (2.58)  and  for  dvaP‘  jd^  -  k“  from  (2.57)  into  (2.59)  leads  to 

^  +  Af+  =  +  A  (v^-  +  -  Q  (V  +  !,“w,  -  f))  ,  (2.60) 

di 

^  +  A1-+  -  k”-^  =  -  Qt"  +  A  (."v*  + 

Q(!/"v  +  3(”'*Wj-1“)),  (2.61) 

which  are  given  in  O’Reilly  and  Turcotte  [52,  eq.  (3.2)].  The  Idnematical  quantities 
associated  with  the  two  motions  are  related  by  (2.54)  and 

r'*’  =  Qr  +  q,  (2.62) 

where  q  represents  a  translation  depending  on  time  alone. 

As  explained  in  O’Reilly  and  Turcotte  [52],  the  fields  f+  and  1"+  cannot,  in  gen¬ 
eral,  be  uniquely  determined  from  (2.54),  (2.62),  (2.60)  and  (2.61),  if  no  invariance 
requirements  are  imposed  on  the  constraint  responses  n,  m"  and  k“.  In  other  words, 
it  is  possible  in  a  constrained  theory  to  have  two  motions  of  the  same  Cosserat  curve 
which  differ  by  a  superposed  rigid  body  motion,  but  whose  constraint  responses,  as¬ 
signed  forces  and  assigned  director  forces  are  not  completely  related.  Furthermore,  the 
boundary  conditions  on  the  respective  n  and  m“  for  the  two  motions  are  not  necessar¬ 
ily  related  either.  In  this  dissertation,  we  impose  the  invariance  requirements  (2.55), 
rather  than  (2.53),  which  is  normally  used. 
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2.5.2  Properly  Invariant  Balance  Laws  for  Approximate 
Theories 

In  this  subsection,  we  give  a  short  account  of  the  properly  invariant  theory  of  O’Reilly 
[50]  as  we  intend  to  apply  it  in  this  dissertation.  First,  recall,  from  the  assump¬ 
tions  (2.2)  and  (2.7),  that  the  tensor  F  is  invertible  and  may  be  uniquely  decomposed 
using  the  polar  decomposition  theorem: 

F  =  RU,  (2.63) 

where  U  is  a  symmetric  positive  definite  tensor  and  R  is  a  proper  orthogonal  ten¬ 
sor.  The  properly  invariant  theory  for  the  approximate  theories  of  a  Cosserat  curve 
require  knowledge  of  the  position  r(^,  t)  and  rotation  tensor  R  =  R(^,t)  of  the  rod 
at  some  material  point  of  the  material  curve  ^  which  is  known  as  the  pivot.  The 
auxiliary  motion  (denoted  by  an  asterisk)  is  constructed  by  removing  the  translation 
and  rotation  at  the  pivot  from  the  total  motion  of  the  rod, 

r*(e,r)  =  R^[r(e,t)-r(6r)]-hR(0,  d:(e,r)  =  R^d„(e,t),  t*  =  t  +  c\  (2.64) 

where  R(0®  refers  to  the  position  of  the  pivot  in  the  reference  configuration  and  c* 
is  a  real-valued  constant.  By  construction,  a  rod  with  small  strain  but  large  motion 
may  have  a  small  (and  possibly  infinitesimal)  auxiliary  motion. 

The  properly  invariant  vectors  n*,  k"*  and  iri“*,  axe  related  to  their  full  motion 
counterparts  by  the  relations 

»•({,(•)  =  R’'nK,i),  rii-K.C)  =  R^m"«.t).  (2.65) 

^We  distinguish  the  referential  position  of  the  Cosserat  curve,  R(^)  (a  vector),  from  the  rotation 
tensor  by  always  showing  its  explicit  dependence  on  ^  alone. 
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The  assigned  forces  f  and  director  forces  1“,  as  noted  previously,  need  not  follow  this 
law  of  transformation,  as  we  do  not  insist  that  the  Lagrange  multipliers  are  invaxiant 
under  superposed  rigid  body  motions.  On  the  other  hand,  it  is  completely  reasonable, 
on  physical  grounds,  to  expect  that  the  direction  of  the  constraint  forces  will  behave 
objectively.  Thus,  recalling  (2.38)  -  (2.40),  we  assume  that’’ 


n  =  =  plR/3“^*,  ih“  =  =  PlRC“^*, 


(2.66) 

where  L  =  1, 2, 3, . . . ,  i2).  The  inertia  terms  associated  with  the  motion  must  be 
accounted  for  in  the  theory,  hence  the  balance  laws  for  the  auxiliary  motion  depend 
on  the  pivot  motion.  However,  it  is  important  to  note  that  the  Lagrange  multipliers 
of  the  auxiliary  motion  are  not  equal  to  pi  as  in  [50]. 

We  will  use  the  auxiliary  motion  to  ensure  that  the  constitutive  laws  are  properly 
invariant.  We  do  so  by  premultiplying  the  balance  laws  for  the  motion  by  RR^  and 
by  using  (2.65)  and  (2.66): 

'dh*  d  (PLV^*) 


ai 


=  A  (v  +  y"w„  -  f ) , 


ft  ^  -pilS”''-  -  k“- 


A  (y"v  +  -  1“)  • 


(2.67) 


(2.68) 


As  in  O’Reilly  [50],  the  constitutive  relations  used  in  (2.67)  and  (2.68)  are  properly 
invariant  under  superposed  rigid  body  motions.  This  invariance  holds  even  if  approx¬ 
imate  constitutive  relations  are  used.  We  note  here  that  O’Reilly’s  developments  are 
based  on  the  work  of  Casey  and  Naghdi  [5]  and  [6], 


'^This  is  equivalent  to  the  invariance  requirement  (p^  =  <p- 
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2.5.3  Modified  Auxiliary  Motion 


Since  the  auxiliary  motion  of  O’Reilly  requires  knowledge  of  the  pivot  rotation  tensor 
R(^,t)  and  displacement  it  is  often  difficult  to  accommodate  it  in  applica¬ 

tions,  including  those  discussed  in  this  dissertation.  In  the  case  of  motions  with  an 
unknown  pivot  rotation  tensor  R(^,t)  (which  arises  particularly  in  flexural  defor¬ 
mations  where  R(|^,t)  is  not  specified  by  the  boundary  conditions),  it  is  helpful  to 
introduce  a  rotation  tensor  S{t)  as  a  known  rotation  such  that  S^S  =  I,  det(S)  =  1 
and  S"*"  =  QS.  When  such  a  rotation  tensor  can  be  found,  a  modified  auxiliary  mo¬ 
tion  can  be  constructed  by  the  same  approach  as  in  [50].  The  difference  between  this 
modified  auxiliary  motion  and  the  usual  auxiliary  motion,  constructed  from  a  known 
pivot  rotation,  is  a  rigid  body  rotation;  and,  therefore,  all  of  the  conclusions  made  by 
O’Reilly  regarding  the  auxiliary  motion  are  equally  valid  for  this  modified  auxihary 
motion.  The  modified  auxiliary  motion  |f,  is  defined  in  a  similar  manner  to 
(2.64):« 

=  S’-(i)  (r((,t)  -  r(f,f))  +s(«),  =  S^{^)d.(^,f).  (2.69) 

where  s{t)  is  a  specified  function  of  time,  t  =  t  +  c  and  c  is  a  constant.  Balances 
of  linear  and  director  momentum  that  use  the  modified  auxiliary  motion  can  be 
constructed  in  a  manner  paralleling  the  development  of  (2.67)  and  (2.68). 


®In  a  personal  communication,  Professor  Casey  informed  Professor  O’Reilly  that  a  related  con¬ 
struction  was  used  in  an  unpublished  draft  of  Casey  and  Naghdi  [5].  The  construction  is  also  related 
to  several  others  in  the  literature,  for  example,  Kane  and  Levinson  [35]  and  Simo  and  Vu-Quoc  [58]. 
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Chapter  3 


Kinematics  of  Approximate  Rod 
Theories 

In  addition  to  the  full  nonlinear  rod  theory,  approximate  rod  theories  can  be  de> 
rived  depending  on  the  magnitudes  of  the  strain  and  rotation  in  the  body.  Casey 
and  Naghdi  [6]  have  developed  related  approximate  theories  for  the  case  of  three  di¬ 
mensional  elasticity.  In  the  case  of  rods,  the  development  is  complicated  by,  among 
others,  the  additional  variable  K^.  Four  approximate  rod  theories  are  developed 
in  the  following  sections.  These  are  infinitesimal  strain  with  infinitesimal  rotation, 
infinitesimal  strain  with  moderate  rotation,  moderate  strain  with  infinitesimal  rota¬ 
tion  and  moderate  strain  with  moderate  rotation.  Before  detaihng  these  theories,  we 
wish  to  make  further  kinematical  developments  of  a  general  nature  beginning  with  an 
equation  for  the  rotation  tensor  R.  The  corresponding  equation  in  three-dimensional 
continuum  mechanics  was  first  derived  by  Shield  [56],  and  its  counterpart  for  Cosserat 
surfaces  was  derived  by  Naghdi  and  Vongsarnpigoon  [49]. 

Again  recalling  the  polar  decomposition  theorem  F  =  RU,  where  U  is  a  symmetric 
positive  definite  tensor  and  R  is  a  proper  orthogonal  tensor,  we  use  this  decomposition 
to  derive  a  differential  equation  for  the  rotation  tensor  R  in  terms  of  U  and  the  other 
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kinematic  variables  defined  in  section  2.1.  We  accomplish  this  by  first  observing  that, 
in  view  of  (2.63),  the  orthogonality  of  R  and  the  symmetry  of  U, 


(3.1) 


Substitution  of  (2.10)  into  the  left-hand-side  of  (3.1)  yields^ 


=  5U-'  [(K,  +  KI  -  2E  „G„  -  2„GIe)  (Ds  ®  D“)  - 
2 


(D“  (2)  D3)  (k„  +  -  2E  oGc  -  2oG^E)  + 


^  (Ds  ®  D“)  -  (D=  ®  Ds) 


d(  +  d( 


dC 


V-\  (3.2) 


This  equation  relates  the  Cosserat  curve  strain-like  deformations  to  the  corresponding 
rotations.  Using  results  developed  by  Hoger  and  Carlson  [31],  the  tensors  U  and  U 
can  be  expressed  in  terms  of  the  tensor  C  =  2E  +  I  and  principal  invariants  of  U : 


Iu  =  trV,  IIu  =  ^[{trVf-trV^],  UIu  =  detlJ.  (3.3) 


Reference  [31]  also  provides  formulae  for  relating  these  invariants  to  those  of  the  tensor 
C.  We  refrain  from  writing  the  lengthy  expressions  here  and  refer  the  reader  to  their 
paper  for  the  results.  In  addition,  the  tensor  did  jd^  can  be  expressed  in  terms  of  U 
(  E  )  and  5E/9^  using  another  residt  of  Hoger  and  Carlson  [32,  eq.  (2.3)]:^ 


dV 

di 


di 


K  J 


{IuIIu-IIIu)niu 
{lullu  -  IIIu)  “ 

Illh  +  [lilllu  +  (luIIu  -  IIIu)  IIu]  .  (3.4) 


^The  result  (3.2)  can  be  verified  by  direct  substitution  of  the  definitions  in  section  2.1). 
*Guo  [34]  concurrently  derived  related  results  in  alternate  forms. 
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Given  E(^,<)  and  Ka(^,t),  (3.2)  can  be  used  to  solve  for  R(^,t)  or,  as  we  intend 
to  show,  as  a  means  of  relating  the  size  of  certain  parameters  of  the  rotation  tensor 
to  these  other  measures.  These  size  relationships  axe  the  basis  of  categorizing  and 
estabhshing  the  approximate  theories  we  will  develop. 

Naghdi  and  Vongsarnpigoon  [49]  (cf.  also  Vongsarnpigoon  [63])  have  shown  that 
the  solutions  R  of  an  equation  that  is  the  equivalent  of  (3.2)  for  a  Cosserat  surface 
may  differ  at  most  by  a  rigid  body  rotation.  The  proof  of  the  corresponding  result  for 
the  case  of  Cosserat  curves  parallels  their  proof  closely,  and  we  do  not  provide  it  here. 
We  record  here  for  future  use  the  well  known  representation  of  R  on  its  eigenbasis: 


R  =  Hs  +  cosO  (g)  /tj  +>2  ®  A*2)  +  sin  0  (S)  A»i  -  Mi  ®  M2) »  (3-5) 


where  Ms  =M3(^)  0  is  the  unit  vector  in  the  direction  of  the  eigenvector  of  R  that  has 
a  unit  eigenvalue,  Ma(^>f)  two  unit  vectors  in  the  plane  perpendicular  to  Ms  such 
that,  for  each  (Mi)/*2>A*s)  forai  a  right-handed  orthonormal  basis  for  5^,  and  6  is 
the  angle  of  rotation.  From  (3.5),  it  follows  that 


jdR 

2 

=  2 

dd 

2 

+  4 

aMs 

03 

<rn 

(1  —  COS0)  . 


(3.6) 


This  result  can  be  estabhshed  by  a  long,  but  standard,  calculation  which  employs  the 
orthonormality  of  fti. 

Detailed  analysis  of  the  approximate  theories  requires  the  expression  of  the  kine¬ 
matic  relations  in  component  form.  In  preparation  for  this,  we  now  introduce  the 
additional  kinematic  variables 


Si  =  d,  -  Di, 


(3.7) 
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and  their  components  as 


6ij  =  6iBj.  (3.8) 

We  now  proceed  to  express  the  kinematic  tensors  E  and  in  terms  of  the  vari¬ 
ables  (3.8).  First,  without  loss  in  generality,  we  assume  that  the  set  {D,}  is  orthonor¬ 
mal:  D‘  =  D,  and  D,  •  Dj  =  Sij.  Then,  substituting  from  (3.7)  into  (2.16)  and 
using  (3.8),  the  tensor  components  Eij  are 

Eij  =  -  -f  6ji  +  SikSjkj  .  (3-9) 


Similarly,  the  non-trivial  components  of  K®  are 

(d.  ■  (d.  •  (d. 


Note  that,  for  initially  straight  rods,  for  which  D,  are  constant,  Eij  are  unchanged 
but  Kai  simplify  to 


Kcgi  — 


dSqj  dSgj  , 


(3.11) 


To  make  the  concepts  of  moderate  and  infinitesimal  theories  precise,  we  define  three 
measures  of  smallness: 


dE 


(3.12) 


3.1  Infinitesimal  Theory 

In  the  infinitesimal  theory,  having  assumed  that  E  is  of  O(eo),  as  cq  — >  0  we  retain 
only  terms  of  this  order  throughout  the  theory.  We  now  introduce  the  concept  of  an 
infinitesimal  rotation  using  a  notion  presented  by  Naghdi  and  Vongsarnpigoon  [49] 
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(cf.  also  Vongsarnpigoon  [63]).  The  rotation  tensor  R,  is  said  to  be  infinitesimal  in 
comparison  to  the  measure  E  if,  for  any  unit  vector  v, 

(R  —  I)  V  =  O(eo)  35  Co  — ^  O'  (3.13) 

Since  Rv  -  v  tahes  its  maximum  value  when  v  is  perpendicular  to  /X3,  we  let  v  = 
in  (3.5),  which  yields 

(R  -  I)  =  (cos  ^  -  1)  /Xj  +  sin  B  /Xj-  (3-14) 

After  considering  the  series  expansions  of  the  sine  and  cosine  functions,  it  is  clear 
that,  for  small  0,  (R  -  I)v  has  the  same  order  as 

0  =  O{eo)  as  eo  ^  0.  (3.15) 


The  approximate  formulae  for  U  when  E  =  O(eo)  as  eo  — ^  0  are 

U«I  +  E  =  I  +  0(€o),  I  -  E  =  I  -  O(eo), 


ly  =  3  +  O(eo),  IIu  =  3  +  O(eo),  IIIu  =  +O(co)  as  eo  — >•  0.  (3.16) 

The  approximations  (3.16)i,2  are  well  known  in  continuum  mechanics  (cf.,  e.g.,  [5]  or 
[50]),  and  the  remaining  approximations  follow  easily.  The  additional  approximation 


dV 


dE 

di 


+  O(eo) 


(3.17) 


is  obtained  from  (3.4)  by  using  (3.16)1,3,4,5.  The  result  (3.17)  agrees  with  the  approx¬ 
imation  of  Naghdi  and  Vongsarnpigoon  [48,  eq.  (4.12)],  which  was  obtained  using  a 
Taylor  expansion. 
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Employing  the  approximations  (3.16)  and  (3.17)  in  (3.2)  yields 


^  [(k«  +  K^)  (Ds  ®  D“)  +  (D“  ®  Da)  (k„  +  K^)  + 


E  (D“  ®  Da)  (K„  +  Kl)  -  (K„  +  K^)  (Da  ®  D“)  e]  +  ^  (Da  ®  D^)  + 
(D“  0  Da)  [e  oG„  +  oG^E  +  ^  (Ka  +  K^)  e]  -  (D"  0  Da)  ^  - 
E  oG„  +  oG^E  +  ^E  (Ka  +  K^)]  (Da  0  D")  +  0{el)  as  €o  ^  0.  (3.18) 


Taking  the  norm  of  both  sides  of  (3.18)  and  using  (3.12)  reveals  that 


aR 


<  4eo  +  2ei  +  2e2. 


(3.19) 


We  now  assume  that  ej  <  Ci^o  ^2  ^  ^260)  where  Ci  and  C2  are  constants.  This 


leads  us  to  conclude  that 


R"^—  =  O(eo)  as  cq  ^  0. 


(3.20) 


Using  this  result  and  (3.6),  it  is  clear  that 


de 


—  =  0{eo)  as  €0  ^  0. 


(3.21) 


Next,  using  standard  inequalities  for  integrals,  one  can  estimate  the  magnitude  of  0 


as 


!««,<) -»({,i)|</-' 


89 


d<i> 


d<j>  = 


89 


di 


L  =  0{eo)L  —  0{eo)  as  cq  ^  0,  (3.22) 


^0 


where  (^0  and  ^  axe  points  along  £  and  L  is  some  integration  length  along  £  which  is 
at  most  the  full  length  of  the  material  curve.  By  employing  the  invariant  theory  of 
O’Reilly  [50],  the  angle  9{^,t)  is  made  to  vanish.  Thus,  we  may  conclude  that 


9{^,t)  =  O(eo)  as  Co  — *•  0, 


(3.23) 
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where  we  have  used  (3.15).  Although  we  do  not  use  (3.23)  in  the  infinitesimal  theory, 
its  counterpart  in  the  moderate  rotation  theory  is  pivotal. 

The  director  displacements  Si  are  now  approximated  by 

Si  =  (F  -  I)D.-  =  (R  -  I)Di  +  O(eo)  as  eo  0,  (3.24) 

where  we  have  used  (2.63)  and,  subsequently,  (3.16)i.  By  setting  v  =  D,-  in  (3.13) 
and  using  (3.24),  we  see  that 

Si  ■  Dj  =  Sij  =  0{€o)  as  eo  ^  0.  (3.25) 

Clearly,  from  our  ajssumptions  following  (3.19),  62  <  C^eo  as  eo  — »•  0,  where  C3  is  a 
constant.  It  follows  from  (3.12)  and  (3.10)  that 

^  =  0(€„).  (3.26) 

Application  of  these  results  to  the  relations  (3.9)  and  (3.10)  yields 

(3.27) 

These  relations  are  recognized  to  be  the  usual  strain-displacement  relations  of  the 
infinitesimal  Cosserat  curve  theory  (cf.  Green,  Naghdi  and  Wenner  [26]). 

Note  that  setting  O{eo)  =  0(ei),  in  this  case,  is  equivalent  to  assuming  that  the 
orders  of  the  kinematic  variables  Sij  do  not  change  upon  taking  their  partial  derivatives 
with  respect  to  This  assumption  is  normally  made  in  the  infinitesimal  theory 
without  a  clear  justification,  but  in  a  nonlinear  theory  it  requires  closer  scrutiny. 
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3.2  Moderate  Rotation  Theory 

We  are  concerned  here  with  a  kinematical  development  in  which  E  is  infinitesimal 
while  R  is  moderate.  In  this  theory  we  neglect  terms  of  0{ej^)  as  cq  0.  The 
rotation  tensor  R  is  said  to  be  moderate  in  comparison  to  E  (recall  that  E  is  O(eo) 
as  eo  — >  0)  if,  for  any  unit  vector  v, 

(R  —  I)  V  =  0(^0 )  2is  Co  — >■  0.  (3.28) 

Following  the  same  reasoning  as  in  the  previous  section,  we  conclude  that 

^  =  0(e|)aseo-^0.  (3.29) 

Again  parallelling  the  approach  of  the  previous  subsection,  we  combine  (3.20)  and 
(3.22)  with  (3.29)  to  get  the  result 

0{eh  =  0{ey)L  =  0{e,).  (3.30) 

We  now  address  the  size  of  the  displacement  components  %  and  their  partial 
derivatives  with  respect  to  ^  for  this  theory.  Recalhng  (3.24)  and  using  (3.28),  we 
find  that 

^.•  =  0(e|),  4  =  0(e|)aseo^O.  (3.31) 

But  this  is  vcdid  only  for  i  ^  j,  since  (3.9)  and  (3.12)i  can  be  used  to  show  that 
8ij  =  O(eo)  as  Co  —>■  0  for  i  =  j.  Next,  as  =  Ka.D’  (gi  =  0(ei)  as  ei  0,  (3.10) 
gives  us  that 

=  0(ei)  =  0(e^ )  as  eo->  0.  (3.32) 
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But  we  assume  that  £2  ^  as  eo  — ^  0,  where  C4  is  a  constant,  so  that  (3.12)3  and 
(3.9)  reveal  that  dSijjd^  =  0{eo)  as  cq  0  when  i  =  j.  Furthermore,  as  8623! 
is  moderate  but  5^23/5^  is  small,  taking  the  derivative  of  (3.9)  with  respect  to  ^ 
shows  that  8632! be  moderate.  The  same  argument  can  be  made  to  show 
that  8S3iI8(  is  moderate. 

In  summary,  using  the  present  definition  of  moderate  rotation  and  certain  assump¬ 
tions  regarding  the  size  of  various  tensors,  we  have  determined  that 

Sii  =  0{€o)  as  €0  -*  0,  (no  sum  on  i)  (3.33) 

Sij  =  O(e^)  as  eo  0.  (i  7^  j)  (3.34) 

=  O(eo)  as  eo  ->  0,  (no  sum  on  i)  (3.35) 

^  =  0(e|)  as  eo  0.  {i  7^  j)  (3.36) 

8^ 

These  results  allow  for  an  easy  simphfication  of  the  strain-displacement  relations  (3.9) 
and  (3.10). 

We  record  here  the  individual  components  of  strain  and  curvature  for  the  initially 
straight  rod  taking  into  account  the  aforementioned  orders  of  magnitude: 


En  =  ^11  +  2  ( 

^12^12  +  ^13^13)  +  ^(^0  )> 

(3.37) 

E22  =  ^22  +  2  ( 

^21^21  +  ^23^23)  +  O(e^), 

(3.38) 

E33  =  ^33  +  —  ( 

^31^31  +  ^32^32)  +  ^(^^)) 

(3.39) 

Ei2  =  E21  =  — 

^^12  +  ^21  +  ^13^23)  +  O(e^), 

(3.40) 

Ei3  =  E31  =  - 

(^13  +  ^31  +  ^12^32)  +  0{eQ ), 

(3.41) 
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E23  =  ^32  =  “  (^23  +  ^32  +  ^21<^3l)  +  ^^(^0  )?  (3.42) 

«„  =  ^  +  ^Ji2  +  ^5.3  +  0{(S),  (3.43) 

di  di  di 

««  =  +  §1S:.  +  +  0(ej),  (3.44) 

di  di  di 

K12  =  +  ^^^23  +  0{€^ ),  (3.45) 

oj  oj 

/C2i  =  ^  +  ^^i3  +  0(e|),  (3.46) 

«i3  =  ^  +  ^^32  +  0(6|),  (3.47) 

dj,  dj 

«23=^  +  ^^3i  +  0(e|),  (3.48) 


as  eo  — 0.  It  is  important  to  note  that  the  length  of  I  poses  restrictions  on  the 
applicability  of  this  theory  through  the  measures  Ci  and  e^. 

3.3  Moderate  Strain  With  Small  Rotation 

Here  we  concern  ourselves  with  a  kinematical  development  in  which  E  is  moderate 
while  R  is  at  most  infinitesimal  (i.e.,  we  let  9  be  of  O(eo)  as  cq  0). 
theory,  we  neglect  terms  of  0{el)  as  cq  0.  Under  these  assumptions,  the  logic  we 
used  for  the  first  two  cases  breaks  down:  we  cannot  use  the  rotation  tensor  equation 
(3.2)  to  estimate  the  order  of  the  displacement  components.  Instead  we  revert  to 
a  geometrical  approach,  recognizing  that  the  shearing  motions  axe  associated  with 
changes  in  the  angles  between  the  directors  whereas  the  rotation  is  associated  with 
the  mean  rotation  of  the  directors.  To  keep  the  rotation  small  while  allowing  the 
strain  to  be  moderate  we  cannot  merely  say  that  some  of  the  components  Sij  are 
moderate  while  others  are  small.  In  fact,  the  remarks  made  above  regarding  the 
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director  motions  reveals  the  shortcomings  of  the  general  approach  we  have  taken.  As 
it  turns  out,  those  shortcomings  are  not  critical  for  the  case  of  moderate  rotation 
with  small  strain,  but  an  improved  approach  is  needed  for  the  case  considered  in  this 
section.  We  detail  such  an  approach  in  the  next  chapter  and  do  not  attempt  to  obtain 
any  strain-displacement  relations  for  this  case  in  the  current  section. 

3.4  Moderate  Strain  and  Rotation 

Again,  we  cannot  use  the  approach  we  employed  in  section  3.2,  but  it  should  be 
apparent  that  all  of  the  components  Sij  will  be  moderate.  In  this  case,  however,  we 
need  to  retain  terms  up  to  O(co)  as  eo  — ^  0.  This  results  in  retaining  the  full  nonlinear 
strain- displacement  relations  (3.9)  and  (3.10). 
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Chapter  4 


Another  Approach  to  Approximate 
Theories 

In  this  approach,  which  parallels  the  work  of  Casey  and  Naghdi  [6],  no  estimate  of 
the  order  of  the  angle  of  rotation  t)  or,  consequently,  the  individual  components 
Sij  is  made.  Rather,  estimates  of  the  orders  of  their  sums  and  differences  are  made. 
This  results  in  a  slightly  more  complicated  expression  for  the  strain-displacement 
relations  in  the  moderate  rotation  theory.  However,  as  will  be  seen  from  the  following 
developments,  the  approach  is  straight-forward.  We  will  also  show  that,  for  the  case 
of  small  strain  accompanied  by  moderate  rotation,  it  can  be  made  to  agree  with 
the  previous  results  to  the  order  of  approximation.  The  resulting  relations,  while 
appearing  significantly  different,  are  expected  to  yield  similar  results  for  all  but  a 
small  class  of  deformations. 

We  begin  by  introducing  the  relative  deformation  measures 

H  =  F  -  I  =  0  D‘  =  e  +  w,  (4.1) 

where 

e  =  i  (H  +  H^)  =  =  i  (J,,  +  Si,)  (Di  ®  D')  .  (4.2) 

w  =  i  (h  -  H^)  =  -w^  =  i  (4  -  Si.)  (d'  ®  D')  .  (4.3) 
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The  tensors  e  and  w,  as  in  the  case  of  three-dimensional  elasticity,  are  associated  with 
the  strain  and  rotation  respectively;  but  this  definition  is  different  from  Naghdi  and 
Vongsaxnpigoon  [49]  at  least  in  form.  In  the  approach  of  this  chapter,  we  detail  four 
approximate  theories  according  to  the  relative  orders  of  these  variables.  For  example, 
the  infinitesimal  theory  corresponds  to  the  case  where  both  e,  w  and  their  parital 
derivatives  with  respect  to  ^  are  small;  infinitesimal  strain  with  moderate  rotation 
corresponds  to  small  e  but  moderate  w  and  5w/5^;  etc.  We  begin  by  defining 


de 

e3{t)  =  sup  II  e  II, 

II 

e5(t)  =  sup  , 

ee{t)  =  sup 

(4.4) 

We  should  point  out  here  that  we  are  assigning  smallness  measures  to  tensor  quantities 
that  are  not  properly  invariant.  This  departs  from  the  approach  of  the  previous 
chapter,  but  we  correct  the  situation  by  the  use  of  the  auxiliary  motion,  which  we 
discuss  further  necir  the  end  of  this  chapter.  In  terms  of  the  variables  e  and  w,  the 
primary  kinematic  tensors  are  given  by 

E  =  e  +  -  +  ew  —  we  —  ,  (4.5) 

2 


K„  =  - 


oG„  =  (2e  -  we  +  +  ew  -  w^)  (8)  + 


de  dw  de  dw  de  5w\  3  .  . 


4.1  Infinitesimal  Theory 

As  stated  previously,  in  this  approach  the  infinitesimal  theory  is  developed  by  assum¬ 
ing  that  all  deformation  measures  are  infinitesimal.  The  appropriate  mathematical 
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statement  of  this  assumption  is 


ce  ^  Cs€5  <  C&ti  <  CtCs  as  C3  0.  (4.7) 

in  (4.7),  C5,  Ce  and  C7  aie  constants.  With  this  assumption,  the  kinematic  measures 
are  given  by 

E  =  e  +  O(e|)ase3-^0,  (4.8) 

K„  =  2e  (g)  ®  +  ^(^3)  as  63  ^  0.  (4.9) 

These  are  again  the  strain-displacement  relations  of  the  infinitesimal  Cosserat  curve 
theory  (c/.  [26]  and  section  3.1).  The  components  of  Ka  given  in  the  previous  chapter 
can  be  recovered  by  substituting  for  e  and  w  and  using  the  identity  D„  •  dD^/d^  = 


4.2  Infinitesimal  Strain  With  Moderate  Rotation 


We  now  assume  that 


65  ^  C'sCe  ^  ^*964  <  ^"1063  as  C3  ^  0, 


(4.10) 


where  Cg,  C9  and  Cio  are  constants.  Then  the  application  of  (4.4)  to  (4.5)  and  (4.6) 
yields 


E  =  e  —  -w^  4-  0(63 )  as  63  — ^  0, 


(4.11) 


K„  =  (2c  -  w^)  ^  +  0(^1)  as  es  -►  0. 


(4.12) 
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The  components  of  deformation  Eij  =  D,-  •  EDy  and  Kai  —  D,-  •  (Kq,D3)  for  an 


initially  straight  Cosserat  curve  using  this  definition  of  moderate  rotation  are 


—  8ki8kj 

—  I  +  0(63  ), 

(4.13) 

88ja\  1 

'  8^  J' 

+  0(e| ), 

(4.14) 

as  ^3  0,  where  we  have  used  (3.9)  and  (3.10)  in  addition  to  (4.11)  and  (4.12). 

Individually,  these  components  are 


En  =  <^11  +  g  (<^31  —  ^13)  +  g  (^21  —  ^12)  +  0{e3 ), 


E22  =  ^22  +  g  (^32  —  ^23)  +  g  (^21  —  ^12)  +  ^^(^3  )> 

E33  =  S33  +  -  ^^32  —  ^23)  +  g  (^31  —  ^13)  +  0(63  ), 


Ei2  =  E21  =  -  ^12  +  ^21  —  ~  (^31  —  ^13)  (^23  —  ^32)!  +  <^(^3  )> 

Ei3  —  E31  —  —  Si3  +  ^31  —  —  ^^21  —  ^12)  (^32  “  ^23^1  +  ^(^3)5 

E23  =  E32  =  -  ^23  +  ^32  —  “  (^21  —  ^12)  (^13  ~  ^3l)  j  +  <^(^3  ), 


(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 

(4.20) 


_  1 

4  V  d( 


S12  —  + 


1  /  5^13  5^31 

4 


^  (<^13  —  ^3i)  +  0{e^ ),  (4.21) 


K22  — 


8822  1  (dhx  88 


8i  '  4:  \8i  8i 


1  /  5^23  ^^32 

4 


5^12  ,  1  ( 


di 

88, 


S23  —  ^32)  +  ^(^3 ))  (‘^•22) 


'31 


di 
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^23 


— ^32) + 0(63 ), 


(4.23) 


«.3  =  ^  +  J  (^  -  ^)  (S3,  -  &)  +  0(4 ).  (4.25) 

(4.26) 

as  63  — >  0. 

We  will  now  show  that,  to  the  order  of  approximation  considered,  this  theory  is 
equivalent  to  that  of  section  3.2.  We  begin  by  recalling  (from  (3.31))  that  the  vectors 
Si  and  their  components  Sij  axe  moderate.  Because  the  Eij  are  infinitesimal,  we  can 
make  the  substitution 

Sji  =  2Eij  -  Sij  =  0(63)  -  6ij  as  €3-^0  (4.27) 

into  (4.15)  -  (4.26)  rendering  them  identical  to  those  of  the  previous  section.  For 
example,  (4.17)  agrees  with  (3.39)  after  letting  €3  <  CnCo,  where  Cu  is  a  constant, 
and  substituting 

S23  =  2E23  —  ^32  =  —  S32  +  0(60),  ^13  =  2E31  —  831  =  —  S31  +  O(eo)  as  60  0  (4-28) 

O  /ij 

into  (4.17)  and  neglecting  terms  of  0{eQ  )  as  cq  0. 

4.3  Moderate  Strain  With  Infinitesimal  Rotation 

Here  we  assume  that 

65  <  Ci2€q  <  C'i3e4  <  Ci4e3  as  €3  0,  (4-29) 
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where  Cn,  C13  and  C14  are  constants.  With  these  assumptions,  the  kinematic  tensors 
axe  given  by 

E  =  e  +  +  0(63)  as  63  0,  (4.30) 

K.  =  (2e  +  e=)  ^  +  ^)  D,  ®  D’  +  0(4)  as  0.  (4.31) 

Since  the  constitutive  coefficients  for  this  theory  are  not  available,  we  do  not  write 
the  component  forms  of  the  strain-displacement  relations. 

4.4  Moderate  Strain  With  Moderate  Rotation 

In  this  theory,  we  assume  the  same  size  comparisons  as  in  the  infinitesimal  theory, 
but  we  retain  terms  of  0(63).  With  this  assumption,  all  terms  of  the  full  nonlinear 
strain-displacement  relations  (4.5)  and  (4.6)  must  be  retained  (the  components  of  the 
relevant  tensors  are  given  by  (3.9)  and  (3.10)). 


4.5  Lack  of  Invariance  of  Approximate  Theories 

As  we  have  stated  previously,  O’ReiUy  [50]  has  shown  that  the  infinitesimal  theory 
is  not  properly  invariant  under  superposed  rigid  body  motions.  We  have  stated  that 
this  is  also  true  for  the  approximate  theories.  We  will  show,  as  an  example,  that  the 
moderate  rotation  theory  is  not  properly  invariant.  Thus  we  have  that 

E'*'  =  e'*'  —  ^  )  as  €3  0.  (4-32) 
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Casey  and  Naghdi  [6]  have  shown  that,  for  rigid  body  motions  superposed  on  a  rigid 
motion  (F  =  1,  E  =  e  =  w  =  0), 

e+  =  -(Q-lf  (Q-I),  w+  =  i(Q-Q^).  (4.33) 

Making  these  substitutions  into  (4.32)  yields 

E+  =  (Q  +  Q’’  -  2l)  +  ^  (21  -  Q"  -  (Q^)")  +  0(e|)  /  0  as  ea  ^  0.  (4.34) 

Clearly  this  theory  is  at  most  invariant  for  superposed  rigid  body  translations  (Q  =  I). 
Thus,  the  moderate  rotation  theory  is  not  properly  invariant  under  arbitrary  super¬ 
posed  rigid  body  motions.  This  situation  can  be  resolved  by  using  the  auxihary 
motion  or  the  modified  auxiliary  motion. 

Omitting  details,  we  now  apply  the  auxiliary  motion  to  this  entire  chapter,  thus 
rendering  the  measures  e  and  w  objective.  In  doing  so,  any  difficulty  with  the  measure 
of  smallness  €3  being  assigned  to  tensors  that  are  not  properly  invariant  is  resolved.  We 
note  further  that,  because  either  the  auxiliary  motion  of  the  modified  auxihary  motion 
fixes  the  rotation  at  the  pivot,  the  solution  R(^,t)  of  (3.2)  is  now  uniquely  specified. 
In  conjunction  with  the  discussion  preceeding  (4.28),  this  places  the  approach  of 
Chapter  3  into  complete  correspondence  with  the  approach  of  this  chapter. 
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Chapter  5 


A  Constrained  Theory  With  Small  Strain 
and  Moderate  Rotation 

We  wish  to  use  the  approximations  E  =  E*  =  and  =  K*  =  O(eo^^^) 

to  simplify  the  equations  governing  the  motion  of  a  Cosserat  curve.  In  order  to 
illuminate  the  desired  simplifications,  we  constrain  the  Cosserat  curve  and  assume  it 
is  initially  straight. 

5.1  A  Particular  Constrained  Theory 

As  we  desire  to  retain  shearing  motions  and  longitudinal  extensions,  we  choose  to 
constrain  the  lateral  extensions.  Thus  we  impose  the  constraints  =  0  and 

=  E22  =  0.  Then,  using  (2.38)  -  (2.40),  we  find  that 

n  =  0,  ic^  =  pidi,  P  =  p2d2,  =  0.  (5.1) 

Substituting  these  expressions  into  (2.33),  and  subsequently  into  the  balance  laws 

(2.67)  and  (2.68)  for  the  auxiliary  motion  gives  us  the  balance  laws  of  the  constrained 
^We  have  returned  here  to  the  asterisk  notation  to  signify  the  auxiliary  motion 
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Cosserat  curve: 


+  Af  =  0,  (5.2) 

R(^^-k^*) +Aq^  =Pidi,  (5.3) 

R  - +  Aq^  =  P2d2.  (5.4) 

By  taking  the  inner  product  of  (5.3)  and  (5.4)  with  d^,  we  now  arrange  the  balance 
laws  in  two  sets:  one  set  having  the  Lagrange  multipliers  absent,^  and  the  other  set 
providing  the  Lagrange  multipliers  in  terms  the  solution  to  the  first  set.  At  the  same 


time,  we  incorporate  the  constitutive  laws  (2.34)  -  (2.36): 


dr  j  .  dr  5d„' 

dfc  +  A- 


'^<3 


'  d 
'  d 


'dr 


djt  + 


+  Af  —  0, 

dr  ^d/jA 

dKlp  d(  ) 


dK* 


dfc  +  Aq^  -  A 


dr  .  dr  ddfl' 

dfc  + 


2k 


dE. 


2k 


dK^U  dC 


or  .  ,  dr  ddu" 

dfc  + 


Oik 


dKl^  d^ 


Pi 


dr  ddfi 

dKlii  d^  ^ 


•d^  -  A 


■  d^  =  0, 
•  d^  =  0, 
.  d^  =  0, 

dr 


P2 


=  ^  A 


dr 


I  df^ 


dfc  +  Aq' 


2k 


_  (dr ddp^ 
d( , 


•  d^- A 


dE*n 

dr 


dE^^ 


(5.5) 

(5.6) 

(5.7) 

(5.8) 

(5.9) 
(5.10) 


In  evaluating  the  partial  derivatives  of  r  tke  balance  laws,  the  partial  derivatives 


axe  first  evaluated  and  then  the  constraints  are  imposed. 

^We  note  here  that  the  complete  decoupling  in  this  constrained  theory  appears  to  be  fortunate 
in  that  the  component  balance  laws  governing  Cosserat  curves  with,  for  example,  shear  constraints 
( El^  =  £^23  =  0)  require  the  substitution  of  the  Lagrange  multipliers  from  the  balance  laws  (this 
can  be  seen  in  the  equations  of  Naghdi  and  Rubin  [47,  §9]  for  a  Bernoulli-Euler  type  rod). 
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5.2  Constitutive  Considerations 


If  we  expand  the  derivatives  of  t/>*  in  (2.34)  -  (2.36)  in  a  Taylor  series  about  the 
origin  and  retain  only  terms  of  O(eg^)  or  larger,  we  would  get  cubic  terms  in  the  cases 
of  moderate  strain  or  rotation.  Although  the  moderate  strain  assumption  results  in 
considerable  simplification  of  the  approximations  for  U(^,t)  and  R(^,t);^  and  the 
associated  kinetical  vectors  n*,k"*  and  the  coefficients  of  the  constitutive  laws 
have  not  been  determined  for  cubic  or  quartic  terms  in  the  free  energy  0*  using 
comparisons  with  the  three-dimensional  theory  of  elasticity,  and  the  constitutive  laws 
are  stiU  nonlinear.  Therefore,  we  will  not  attempt  to  deal  with  the  constitutive  laws 
for  these  cases,  as  there  is  no  apparent  advantage  to  the  approximate  theory  over  the 
full  nonlinear  theory. 

Instead,  we  focus  on  a  specific  free  energy  which  has  appeared  previously  in  the 
literature.  Our  free  energy  will  not  be  general  enough  to  encompass  all  cases  of 
moderate  rotation,  but  the  precise  values  of  its  coefficients  are  known.  The  free 
energy  of  the  rod  that  we  use  is  invariant  under  the  transformations 

d.-^±d.,  D.-^±D.-,  (5.11) 

The  free  energy  !/>*  we  use  is  (adapted  from  [18]  and  [26]):'* 

2A^*  =  4  -|-  OC2E22  +  Q:3-£'33  +  ^5^23  +  <^6^13  -f  Q'7.E*j£?22  d"  ®8.£'il-£'33+ 

a9f^22-^33)  +  *^4  iEl2  +  ■£'2l)^  +  <^10«n  +  «11«22  +  ^121^12  +  <^13'«21  + 

^14^12^21  ^15^23  “I"  ^16^13  +  ^17^11^22  ^0  (5.12) 

^See  [6,  eq.  (3.13)]  for  some  of  these  approximations. 

'^The  notation  ai  —  air  corresponds  to  the  notation  ki  —  kn  in  [18],  [22]  and  [26]. 
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Some  of  the  constitutive  coefficients  for  homogeneous  isotropic  materials  deter¬ 
mined  by  Green,  Laws  and  Naghdi  [19],  Green,  Naghdi  and  Wenner  [26,  §9]  and 
Green  and  Naghdi  [22,  §10]  are  listed  below; 


=  «2  =  «3  =  777 — =  ae  =  EAf{u), 


07  =  Og  =  OC9 


2(1  -  2u) 

fivA 


,  oio  =  fih,  otn  = 


il-2uy 

ais  =  EIi^  ais  ~  El2^  air  =  0, 


(5.13) 

(5.14) 

(5.15) 


where  A  is  the  rod  cross-sectional  area,  7i  and  I2  are  the  referential  area  moments  of 
inertia  about  the  lines  =  0  and  6^  =  0  respectively,  /x  =  jE/2(1  -|-  i/)  is  the  sheax 
modulus,  E  is  Young’s  modulus,  u  is  Poisson’s  ratio  and  f{v)  is  a  function  of  v  which 
has  several  alternative  forms  in  the  literature.® 

Having  completed  the  constrained  theory,  we  should  point  out  that  this  model 
is  stiffer  in  extension  than  a  rod  that  is  unconstrained.  Thus,  such  a  model  will 
underpredict  the  extensional  response  of  a  rod  that  allows  lateral  extension  (this  also 
results  in  higher  natural  frequencies  for  extensional  vibration).  We  could  adjust  the 
value  of  ora  (say  to  equal  EAjA)  to  correct  this  problem,  but  we  let  this  issue  stand 
here. 


5.3  Balance  Laws  for  Constrained  Straight  Rods 

In  this  section  we  assume  the  rod  is  initially  straight  and  incorporate  the  free  energy 

expression  (5.12)  as  well  as  the  strain-displacement  relations  (3.37)  -  (3.48)  to  reduce 

®[26,  eq.  (9.29)]  cind  [22,  eq.  (10.29)]  have  slightly  different  values  for  f{v)  which  were  obtjuned 
by  different  procedures. 
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the  balance  laws  to  a  set  of  seven  equations  in  the  seven  unknown  displacement 
variables.®  To  do  this,  we  first  note  that 


q”  =  0(e;),  !  =  0(6;),  (6.16) 

d“  =  D|  +  0(4’),  A,:  =  ^  +  0(e;)  as  4  -*  0.  (5.17) 

We  now  change  to  the  tilde  notation,  which  refers  vectors  to  the  corotational  basis 
|RD,}.  On  this  basis,  the  balance  laws  (5.5)  -  (5.8)  become,  ^ 

^  (^13  +  ^31  +  ^12^32)  +  «i5  +  ^~f  =0,  (5.18) 

^  (^23  +  ^32  +  ^2i4i)  +  “16  9^)  ^  ~  (5.19) 

2  2\ 

^  ^40(3^^33  +  ^  (^31^31  +  ^32^32)^  +  ^16  (~^)  j  ~ 


d  f  ( d8i2  5^13  ?  \  .  1  ( ^^21  ,  ^^23  c  ^  ^  I  d6z2 

^  (  ai2  ( -^  +  -^<^23j  +  I  —  +  —^13  1  1  +  o;i6— —  + 


2””  \di  '  dt  ~^yj  '  d( 

—  a4  ^^12  +  ^21  +  ^13^23)  =  O5  (5-21) 


d  f  ( 8821  ,  3^23 1  \  , 

a?  (it  +  IT  “j 


5““( 


^<^12  .  ^^13  r  \\  .  ^^23^^31  , 

+  ~T:ir^2z  11+  :Tr'  + 


di  di 


di  di 


23 


d  (  ( d8\^  ^  5^12  J  W 

8?  l““  VW  ^  'W  VJ  ■  'W 


1  +  ^13^23)  —  0, 

(5. 

f  1  ^ 

d~8i2\ 

-  [a,2  -  -0(14  j 

dU 

+ 

—  06  (^13  +  ^31  +  ^12^32)  =  9,  (5.23) 


®We  do  not  include  the  equations  for  the  Lagrange  multipliers. 

■^The  kinematic  variables  are  the  same  in  either  basis,  but  the  asterisk  notation  is  not  valid  for 
the  components  of  f  or  1"  in  the  constrained  theory,  so,  for  the  sake  of  consistency,  we  use  the  tilde 
notation  for  cdl  variables. 
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—  Q!5  (^23  +  ^32  +  ^21^3l)  =  0>  (^•24) 


where  we  have  omitted  terms  of  0(e^^^)  as  eo  0-*  Note  that  eo  =  ej.  These 
equations  must  be  supplemented  by  the  displacement  equations  (from  (2.15)) 


-  du2  ~  du3  I 

—  031»  —  032>  —  O33, 

OK 


(5.25) 


and  the  two  constraints: 


Su  —  (6x2^12  +  ^13^13)  )  ^22  —  (^21^21  +  ^23^23)  • 


(5.26) 


With  the  constrained  equations  in  this  form,  one  can  see  that  the  static  coupling 
in  the  moderate  rotation  Cosserat  curve  theory  is  minimal.  In  fact,  in  the  absence  of 
flexure  (^13  =  ^31  =  ^23  =  ^32  =  0),  both  the  extensional  and  torsional  equations  are 
completely  hnear.  The  only  couphng  is  between  flexure  and  extension,  and  this  is  only 
a  one-way  couphng  -  the  flexural  equations  being  completely  hneax,  but  the  flexural 
response  providing  a  forcing  input  to  the  extensional  equations.  In  Chapter  7  we 
provide  examples  to  illustrate  these  statements  as  well  as  to  show  that  the  assumptions 
made  in  the  theory  are  vahd  for  some  rod  geometries  and  deformations. 


5.4  Infinitesimal  Theory  of  Cosserat  Curves 

In  this  section  we  record  the  component  equations  for  the  unconstrained  infinitesimal 

Cosserat  curve  theory  for  future  use.  These  equations  are  well  established  in  the 

®In  general,  the  components  /*  and  /"*  may  be  difficult  to  determine,  but,  in  this  dissertation, 
we  restrict  our  attention  to  problems  where  either  f  and  1"  vanish  identically  or  R  =  I. 
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literature,  having  appeared  previously  in  [18],  [20],  [22]  and  [26].  The  assumptions 
that  the  rod  is  initiciUy  straight  and  that  both  E*  and  K*  are  small  lead  to  the  hnear 
balance  and  constitutive  laws  of  Green,  Naghdi  and  Wenner  [26]. 

Due  to  the  assumptions  made  above,  which  result  in  leaving  the  kinetical  vectors 
n*,  m"*  and  k“*  of  O(eo)  as  0  and  A^.  =  0,  the  balance  laws  (2.24)  -  (2.27) 
reduce  to® 

— —  -f-  A/  =  0{el^)  as  ej  0,  (5.27) 

+  xr  =  0{6f)  as  eS  0,  (5.28) 

_  p.A  ^  0^  ^A,3  _  ^*A  ^  0(^*2)  as  ej  ^  0.  (5.29) 

The  constitutive  laws  can  be  obtained  by  the  application  of  (2.34)  and  (2.36) 
to  (5.12),  after  neglecting  terms  of  O(eo^)  as  ej  ^  0.  We  list  here,  in  component 
form,  the  linear  constitutive  laws: 


=  2aeE;3,  n*^  = 

=  4a3£?33  +  2a^Eli  +  2a9£?22? 

(5.30) 

=  Q;io«n  +  ^'^17«22) 

0*0  *  ,  1  « 

m  =  aiiAC22  + 

(5.31) 

m}*^  =  Q:i2Kt2  + 

2*1  *  1  ^  * 

m  =  ai3/C2i  + 

(5.32) 

=  crisKja, 

(5.33) 

=  4.CX1E11  +  20(7^22  +  2oisE'^^,  k‘ 

=  40C2E22  +  2a7Eli  +  209^33, 

(5.34) 

=  404^12. 

(5.35) 

The  hnear  equations  of  motion  for  an  initially  straight,  uniform  Cosserat  curve 

are  obtained  by  combining  the  balance  laws  (5.27)  -  (5.29),  constitutive  laws  (5.30) 

^Note  here  that  because  n  =  Rji*,n'  =  n  •  d‘,  n’*  =  n*  •  d**  and  we  have  an  unconstrained 
theory,  it  follows  that  n*'  =  n‘.  Similarly,  m“*'  =  m“’  and  fc"*’  =  fc"*. 
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-  (5.35)  and  the  infinitesimal  strain-displacement  relations  E*- 


<.•  =  d6*Jdi: 


Q!6^  (^13  +  ^3l)  +  ^  (/^  “  “l)  =  0,  (5.36) 

(^23  +  ^32)  +  A  —  U2)  =  0,  (5.37) 

2q!8-^^  +  2a9-^^  -f  4a3-^^  -f  A  =  0,  (5.38) 

+  -ai4  ^^2^  +  A  -  «4  (^12  +  ^21)  =  (5.39) 

ai3^  +  +  A  -  a4  (^1*2  +  ^2*1)  =  0,  (5.40) 

ai6^  +  A  -  cce  (^1*3  +  4*1)  =  0,  (5-41) 

ai5^  +  A  (P  -  -  «5  (^2*3  +  ^3*2)  =  0.  (5.42) 


Qfio  +  2^17  ^^2^  +  A  j  —  40:1511  —  207622  —  20^633  =  0,  (5.43) 

ail  ^  +  ^"17^  +  A  -  ^""^22)  -  4a25;2  -  2a75ri  -  2a953*3  =  0,  (5.44) 

where  we  have  also  assumed  that  the  coefficients  y^=  y^=  y^^  =  y^^  =  0. 
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Chapter  6 


Free  Linear  Vibration  of  Cosserat  Curves 

In  this  chapter,  we  will  present  examples  from  the  unconstrained  infinitesimal  theory 
of  a  Cosserat  curve.  We  include  the  free  flexural  and  extensional  vibration  of  initially 
straight,  uniform,  isotropic,  rods.  These  results  are  seen  to  be  similar  to  results  from 
other  rod  theories,  but  we  provide  some  new  results  in  the  case  of  flexure  and  more 
general  results  in  the  case  of  extension.  Previous  research  on  the  linear  vibration  of 
Cosserat  curves  include  the  calculation  of  wave  speeds  in  straight  infinite  Cosserat 
curves  by  Green,  Laws  and  Naghdi  [19],  shock  wave  propagation  in  Cosserat  curves  by 
Cohen  and  Whitman  [7],  wave  speeds  in  initially  curved  directed  curves  by  Eason  [14] 
and  extensive  analysis  of  the  flexural  vibration  of  a  Timoshenko  beam  by,  for  example, 
Dolph  [11],  Downs  [12]  and  Huang  [33].  The  latter  vibrations  are  equivalent  to  those 
of  the  infinitesimal  theory  of  flexure  in  directed  curves. 

6.1  Flexural  Vibrations 

In  this  section,  we  address  the  general  problem  of  the  free  flexural  vibration  of  a 
Cosserat  curve  using  the  infinitesimal  theory.  As  a  specific  example,  we  take  a  circular 
cantilever  rod  and  determine  the  modes  of  vibration  in  the  6'^  —  ^  plane.  To  determine 
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these  modes,  we  calculate  the  natural  frequencies  in  which  the  rod  can  vibrate  and 
their  associated  mode  shapes.  By  free  vibration  we  mean  the  response  to  initial 
displacements  and  velocities  in  the  absence  of  applied  loads.  In  simpler  models  of 
beam  vibration,  all  of  the  modes  typically  participate  in  the  free  response;  however, 
for  this  model  (which  may  be  placed  in  correspondence  with  the  Timoshenko  beam 
theory),  there  are  some  modes  that  exist  only  for  special  geometries.  We  consider  all 
possible  modes,  including  one  mode  that  has  been  ignored  in  the  literature  and  was 
reported  recently  by  O’Reilly  and  Turcotte  [51].  For  the  Cosserat  curve  rod  model, 
the  mode  shapes  consist  of  the  curve  displacements  and  director  displacements.  In  the 
case  of  planar  flexure,  there  is  only  the  displacement  U2  and  the  director  displacement 
^23-  These  are  shown  in  Figure  6.1. 


Figure  6.1:  Some  of  the  kinematical  variables  associated  with  the  flexure  of  a 
Cosserat  curve. 

To  ensure  that  the  theory  is  properly  invariant  under  superposed  rigid  body  mo¬ 
tions,  we  construct  the  modified  auxiliary  motion  corresponding  to  the  motion  of  the 
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rod.  Since  the  rod  has  a  fixed  end,  we  choose  this  end  as  the  pivot  =  0).  Then  the 
choice  S  =  I,  ^=0,  c  =  0  and  s{t)  =  0  in  (2.69),  which  represents  the  rotation  of  the 
fixed  object  to  which  the  rod  is  attached.^  We  note  that  the  original  motion  is  identi¬ 
cal  to  the  modified  auxiliary  motion  in  this  case  except  when  rigid  body  motions  are 
superposed  (where  S+  =  Q,  and  the  modified  auxiliary  motion  remains  unaltered  but 
the  original  motion  does  not).  Because  of  the  trivial  nature  of  the  modified  auxiliary 
motion  in  this  example,  we  do  not  use  the  asterisk  or  tilde  notation;  however,  all 
of  the  results  presented  axe  understood  to  be  associated  with  the  modified  auxiliary 
motion. 

The  Cosserat  curve’s  flexural  response  is  governed  by  the  homogeneous  equations 

05^  ^^23  +  -  •'“2  =  0,  (6.1) 

a.3^  -  C.3  +  ^)  =  -'!'“523.  (6.2) 

As  remarked  previously  by  Green,  Laws  and  Naghdi  [18],  these  equations  are  identical 
in  form  to  those  of  the  Timoshenko  beam.  The  solutions  to  these  equations  must 
satisfy  the  appropriate  boundary  conditions  as  well  as  the  initial  conditions 

M(,o)  =  hiO,  fe(S,0)  =  «2(«,0)  =  L((,0)  =  S2(fl.  (6.3) 

Adapting  the  approach  of  Huang  [33]  for  the  Timoshenko  beam,  we  rearrange  these 
coupled  equations  into  two  separate  equations  in  U2  and  ^23-  By  taking  appropriate 

partial  derivatives  and  solving  the  first  equation  for  ^23  to  substitute  into  the  second 

^For  other  boundary  conditions,  such  as  pinned-free  or  free-free,  finding  a  known  rotation  tensor 
S  requires  more  construction  than  we  provide  here. 
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equation  we  obtciin 


015 


d^U2 


-  A  (j,^^  + 


d‘*U2  ,  ^U2 

as  dt^'^  dt^  ’ 


(6.4) 


where  we  have  switched  to  the  partial  derivative  notation  for  material  time  derivatives. 
Similarly,  by  solving  (6.2)  for  5u2/5^,  taking  appropriate  partial  derivatives  of  the 
result  and  substituting  them  into  (6.1)  we  get 


a"^23 


a^2a<2  +  cts  dt*  dt^  ' 


(6.5) 


We  solve  these  equations  by  first  separating  the  variables  into  functions  of  ^  and  a 
harmonic  function  of  time: 


«2  =  U2{i)  sin(a;t  -  <f>),  S23  =  A23(0  sin(a;<  -  (f)).  (6.6) 


Substituting  this  assumption  into  (6.4)  and  (6.5)  and  canceling  the  common  time 
dependent  factor  we  obtain 

d^U2  oi5A  d^U2  .  ,  2  /'w2Ay22 


<^15' 


+  (jj^X 


{y 


+ 


'  — )^  +  Au;^ 


as/  de 


05 


-  1)  £^2  =  0, 


(6.7) 


with  an  identical  equation  for  A23. 

Following  the  standard  procedure  for  solving  equations  of  the  type  (6.7),  we  as¬ 
sume  solutions  of  the  form  U2  =  Cic^  and  A23  =  ^26^^.  We  substitute  these  as¬ 
sumptions  into  (6.7)  and  cancel  the  common  exponential  term  to  get  a  characteristic 
equation  (which  is  obviously  the  same  for  U2  and  ^23): 

[y'‘^  +  ^)  +  -  1)  =  0.  (6.8) 
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Using  the  quadratic  formula  we  solve  for  the  roots  of  (6.8): 


_ 
P  = 


f 


1 


2 

Au;^ 


,ai5  Os 

ai5  «5 


)± 

-Xio^ 

fXy^^u^ 

1 

4 

V«15 

OC5J 

\  Q!15«5 

ais/ 

4  VttlS  Ofs/  «15 


The  roots  axe  clearly  always  real,  but  fall  into  three  categories: 
Case  I 

Xy'^W  <0:5,  Pi  <  0,  pI>  0, 


(6.9) 


(6.10) 


Case  II 

Xy^^Lo^  =  as,  pI  <0,  pI  =  0,  (6.11) 

Case  III 

Ay^^w^  >0:5,  Pi  <  0,  pI<  0,  (6.12) 

where  pi  is  associated  with  the  minus  sign  in  (6.9)  amd  p2  with  the  plus  sign  in  (6.9). 

For  the  moment,  we  focus  on  case  I.  We  now  let  p,  correspond  to  the  magnitude 
of  the  above  roots  p,  and  apply  Euler’s  formula  for  complex  exponentials  to  recast 

the  assumed  solutions  into  sines,  cosines  and  their  hyperboHc  counterparts.  For  case 

I,  the  solution  is 


U2  =  A  cos  Pi ^  +  Bsinpi^  +  C  coshp2^  +  i?sinhp2^,  (6.13) 

A23  =  Ecospi^  +  Fsinpi^  +  Gcoshp2if  +  ifsinhp2^.  (6-14) 

The  eight  coefficients  introduced  above  axe  not  independent  but  «ire  rather  related 
through  the  original  coupled  equations  of  motion.  Either  equation  can  be  used  to 
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find  the  relationships.  Using  (6.1),  the  relationships  are  found  to  be 


The  four  remaining  coefficients  are  subject  to  four  boundary  conditions.  Application 
of  all  four  conditions  yields  the  characteristic  equation  for  the  natural  frequencies 
of  flexure  of  the  rod  (this  can  be  done  either  by  eliminating  the  coefficients  A  —  D 
or  by  forming  a  homogeneous  matrix  equation  of  the  boundary  conditions  whose 
determinant  must  vanish).  Once  the  natural  frequencies  are  estabhshed,  the  eigen¬ 
functions  can  be  determined  by  application  of  any  three  of  the  boundary  conditions. 
The  free  response  of  a  rod  with  the  boundary  conditions  usually  considered  can  be 
represented  as  an  infinite  sum  in  the  products  of  the  eigenfunctions  and  the  time 
dependent  functions  in  (6.6):^ 


oo  _ 

“2(^,0  =  n  t^2„(0sin(w„t  -  <^„),  S23{(,t)  =  Y,  ‘^23„(0  siii(u;„i  -  (6.16) 

ri=l  n=l 

Finally,  the  coefficient  of  each  eigenfunction  in  the  sum  and  the  corresponding  phase 
angle  are  found  by  application  of  the  initial  conditions  (6.3).  In  general,  the  ini¬ 
tial  conditions  must  be  transformed  to  generalized  initial  conditions  (in  the  modal 
coordinates)  before  the  coefficients  can  be  determined. 

Turning  now  to  a  specific  example,  the  boundary  conditions  for  the  cantilever  rod 


are 


^2(0)  =  0,  A23(0)  =  0,  ^A23(.£')  +  (■^)  ~ 


^Dolph  [11]  has  shown  that  the  problem  is  self-adjoint,  and  therefore  the  eigenfunctions  form  an 
orthogonal  basis  on  the  response  space. 
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We  apply  all  of  the  boundary  conditions  to  eliminate  the  coefficients  in  (6.13)  resulting 
in  the  frequency  equation 


—  sinh P2L  -  sin piL]  ( sin piL  +  —  sinh P2L 

P2  /V  Pi 

^cospiZ  —  —  coshp2l^^  (cospiL  —  7coshp2^) »  (6.18) 


where 

=  -  QsPi 

Xuj^+aspl' 

Once  the  natural  frequencies  are  obtained  from  this  equation,  any  three  of  the  bound¬ 
ary  conditions  can  be  used  to  solve  for  the  mode  shapes: 


r  ^  sinh  P2L  —  sin  pi  L 

U2„=A„  cospi^  -  coshp2^  +  — - r — T - F 

7  cosh  P2L  —  cos  Pi  L 


sinpi^  +  7  —  sinhp2^ 
Pi 


(6.19) 


where  the  coefficient  A„  for  each  mode  is  determined  by  the  initial  conditions  (after 
reuniting  these  functions  with  the  time  dependent  part),  and  the  p,-  are  functions 
of  the  mode  natural  frequency  according  to  (6.9).  This  result  agrees  with  that  of 
Huang  [33].  Similar  eigenfunctions  can  be  written  for  A23  by  applying  the  factors 
in  (6.15)  to  the  coefficients  in  the  above  eigenfunctions,  but  we  do  not  record  these 
functions  here. 

We  have  yet  to  examine  the  results  for  ca^es  II  and  III.  Case  III  is  a  simple 
modification  of  case  I,  since  the  only  difference  is  the  sign  of  p\.  We  simply  replace 
every  occurrence  of  p2  by  ip2 ,  where  i  represents  \/— 1.  This  eliminates  the  hyperboHc 
functions  replacing  them  with  additional  trigonometric  functions;  but  the  procedure 
for  determining  the  modes  is  identical.  Therefore,  we  now  proceed  directly  to  case  II. 

In  case  II,  there  is  only  one  non-zero  root  p^,  so  that  the  assumed  solution  should 
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take  the  form 


U2  =  A'  cospi(  +  B'sinpi^  4-  C'^  +  D', 


(6.20) 


A23  =  E'  cos  Pi ^  +  F' sin  pi(  +  G'(  +  H'. 


(6.21) 


As  with  case  I,  the  above  coefficients  are  related  through  the  balance  laws.  Using  (6.1) 

again,  the  relationships  for  this  case  are  found  to  be 

E'  _  Qs  F'  0:5  ^  _  J_  r*' —  0 

B'  Plan'  A'  Plan'  D'  j/22’ 

with 

/  ^.  1  \  2 

(6.23) 


Application  of  the  four  boundary  conditions  in  this  case  leads  to  a  condition  on  the 
geometry /material  of  the  rod  (the  frequency  of  vibration  has  already  been  determined 
in  terms  of  the  geometry  and  material  according  to  (6.11)).  The  eigenfunctions  can 
then  be  determined  from  some  of  the  boundary  conditions.  To  our  knowledge,  only 
simply- supported  beams  have  previously  been  considered  by  TraiU-Nash  and  Col¬ 
lar  [61]  and  Downs  [12],  and  their  analysis  is  incomplete.  Applying  all  four  cantilever 
boundary  conditions  (6.17)  simultaneously  leads  to  the  condition 

(6.24) 


,  PlL  .  j-  ,  ^  f  \  r  A 

1 - r-smpiL  -f-  -  (  — - h  1  cospiL  =  0, 


2  "  '  2  V  ai5  '  y'^'^as^ 

which  depends  on  the  rod  geometry  and  Poisson’s  ratio  but  not  on  the  shear  modulus 

of  the  material.  The  corresponding  eigenfunctions  are 


and 


t^2(0  =  A' {cos pi^  -  7'sinpi^  -  1) , 


A23(0  =  (-^  [sinpi^  -f  7'  (cospi^  -  1)]  -I-  , 

[Pian  y  ) 


(6.25) 


(6.26) 
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where 


,  _  Qi5  +  y^^Q:5  cos  Pi  L 
^  as  sin  Pi  L 

The  case  II  results  for  the  simply-supported  rod  deserve  special  mention  as  they 
are  different  from  all  other  boundary  conditions  including  the  cantilever  case  given 
above.  The  boundary  conditions  for  a  simply-supported  rod  in  flexure  are 

i7,(0)  =  0,  U,{L)  =  0,  ^(0)  =  0,  ^(£)=0-  (8-27) 

Application  of  the  first  and  third  of  these  conditions  leads  to  A  =  Z?  =  0.  The  other 
two  conditions  lead  to  two  possibilities,  which  are  B  =  0  oi  piL  =  mr  (n  =  1, 2, 3, ...). 
If  we  choose  B  =  0,  this  mode  exists  for  all  materials  and  geometries  and  has  the 
simple  shape 

UiiO  =  0,  A23  =  H.  (6.28) 

This  solution  was  first  obtained  by  TraiU-Nash  and  Collar  [61]  and  was  subsequently 
rediscovered  by  Downs  [12].  The  other  mode,  which  occurs  at  the  same  frequency 
and  has  not  been  reported  on  previously  to  our  knowledge,  will  exist  if,  and  only  if, 

P^L=(^  +  ^Y  L  =  nr  (n  =  1,2, 3, ...).  (6.29) 

This  condition  places  restrictions  on  the  material/geometry  combination  sinailar  to 
the  condition  on  the  cantilever  rod.  The  corresponding  eigenfunctions  for  this  mode 
(when  it  exists)  are 

U2{0  =  -S'sinpi^,  A23(0  =  ^  -Pi  j  cospi^  +  H.  (6.30) 

\P^  y  j 

Further  details  of  these  results,  in  the  context  of  Timoshenko  beam  theory,  can  be 
found  in  O’ReiUy  and  Turcotte  [51]. 
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As  an  illustration,  we  determine  the  cantilever  rod  modes  for  ca^es  I  and  III 
by  fixing  the  material  and  geometry  (we  choose  a  circular  steel  rod  with  v  =  0.3, 
R=0,lm  and  L=1.0m)  and  iterate  the  frequency  equation  (6.18)  to  find  the  natural 
frequencies.  We  substitute  these  frequencies  into  the  eigenfunctions  and  plot  the 
resulting  functions  in  Figure  6.2,  where  we  also  plot  the  corresponding  eigenfunctions 
for  A23.  One  should  keep  in  mind  that  "the  amplitudes  of  the  functions  U2{i)  are 
arbitrary  until  the  initial  conditions  are  invoked,  but  the  relative  amplitude  of  A23  to 
U2  is  fixed  by  the  analysis. 

The  first  eight  modes  are  from  case  I,  and  the  remainder  are  from  case  III.  It  is 
interesting  to  note  that  case  I  modes  have  an  additional  node  for  each  new  mode, 
but  this  is  not  true  for  case  III.  The  associated  natural  frequencies  are  plotted  in 
Figure  6.3,  where  we  have  also  in^cated  the  dividing  line  between  cases  I  and  III 
(i.e.,  the  case  II  frequency). 

We  also  note  that  the  case  II  frequency  does  satisfy  the  frequency  equations  for 
both  cases  I  and  III,  but  is  not  a  mode  of  this  rod  since  the  material/geometry 
combination  used  to  compute  these  modes  does  not  satisfy  the  condition  for  case  II 
vibration  (6.24).  To  satisfy  this  condition  for  plotting  purposes,  we  keep  the  length  at 
one  meter  and  Poisson’s  ratio  at  0.3  and  adjust  the  radius  to  0.0909m.  The  resulting 
mode  shapes  are  plotted  in  Figure  6.4.  Other  possible  combinations  of  length  and 
radius  that  satisfy  (6.24)  are  plotted  in  Figure  6.5. 
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Figure  6.2:  The  first  twelve  flexural  eigenfunctions  of  a  cantilevered  (fixed-free) 
Cosserat  curve  with  /i  =  82.7GPa,  =  7500Kg/m^,  v  —  0.30,  L  =  1.0m  and 

R  =  0.1m.  K'7 


Figure  6.3:  The  first  19  flexural  eigenfrequencies  of  a  cantilevered  Cosserat  curve 
with  n  =  82.7GPa,  p*^  =  TSOOKg/m^,  v  =  0.30,  L  =  1.0m  and  =  0.1m. 


6.2  Extensional  Vibrations 

The  Hneax  equations  for  the  extensional  vibration  of  a  Cosserat  curve  are  unique 
in  their  inclusion  of  transverse  extensional  vibration  for  rods  that  axe  not  circular 
(Mindlin  and  Hermann  [44]  did  include  transverse  inertial  effects  for  circular  rods, 
but  no  theory  other  than  the  Cosserat  theory  of  Green  and  Naghdi  includes  trans¬ 
verse  extensional  vibrations  for  non-circular  rods).  Other  rod  theories  that  consider 
transverse  effects  generally  have  only  one  direction  of  freedom,  requiring  the  trans¬ 
verse  response  to  foUow  the  longitudinal  response  according  to  the  Poisson  effect  (for 
example.  Love’s  correction  [41,  §278]). 
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Figure  6.4:  The  case  II  mode  shape  for  a  cantilevered  Cosserat  curve  with  fi 
82.7GPa,  Pq  =  7500Kg/m^,  t/  =  0.30,  L  =  1.0m  and  R  =  0.0909m. 
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Figure  6.5:  Case  II  length  and  radius  combinations  satisfying  the  existence  condition 
(6,24)  when  =  82.7GPa,  Pq  =  7500Kg/m^  and  i/  =  0.30. 

To  keep  the  theory  properly  invariant  under  superposed  rigid  body  motions,  we 
construct  the  auxiliary  motion.  Since  the  rod  has  a  fixed  end,  we  choose  this  end 
as  the  pivot  (^  =  0).  Then  R  =  I  and  the  original  motion  is  again  identical  to  the 
auxiliary  motion  except  when  rigid  body  motions  are  superposed.  Because  of  the 
trivial  nature  of  the  auxihary  motion  we  suppress  the  asterisk  notation,  but  we  claim 
that  all  results  are  properly  invariant  under  superposed  rigid  body  motions. 

The  equations  for  the  hnear  free  extensional  vibration  are  (assuming  y^=y^  — 

5^11  0622  0633  ...  ,  . 

2q!8"'t^  +  2a9-^— +  403-^-  Aria  =  0,  (6.31) 

di 
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(6.32) 


d^6\i  1  0^622  \ii?  A  K  ox  ox  n 

aio  Q.9  +  n^i7-^ - A2/"<5ii  —  4ai(5ii  —  207622  —  203^33  —  0? 

^0^17  —  4^02622  —  2Q!7^11  —  209^33  =  0.  (6.33) 

We  intend  to  use  a  rectangular  section  in  our  example  so  that  0!i7  =  0  (see  (10.42)  of 
[22]).  As  usual,  we  proceed  by  separating  variables  and  assume  that 

U3{(,t)  =  U3{()sm{ojt  -  <^),  ^ii(^,^)  =  Aii(^)sin(a;t  -  (f>), 

S22(C  t)  =  A22(0  sm(wt  -  <^).  (6.34) 

We  also  need  to  incorporate  (5.25)  to  reduce  the  problem  to  three  equations  in  three 
unknowns.  Incorporating  these  eissumptions  and  putting  the  equations  in  matrix  form 
yields 
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Now  since  the  center  matrix  in  the  above  system  cannot  be  made  proportional  to 
the  other  two,  we  multiply  through  by  the  inverse  of  the  first  matrix  and  transform 
the  second  order  system  into  a  larger  first  order  system  using  the  usual  state-space 
approach  (in  our  case,  however,  the  independent  variable  is  spatial  rather  than  tem- 
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(6.36) 


or 


{x')  =  [A]  {x) . 


Ttere  axe  maxiy  methods  in  the  controls  and  linear  algebra  literature  for  solving 
(6.36),  but  our  intent  is  simplicity  rather  than  efficiency.  We  take  the  straight-forward 
approach  of  assuming  a  solution  of  the  form  {x}  =  {C}  The  general  solution  that 
follows  from  this  approach  may  be  written  as 
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(6.37) 


i=l 


where  the  pi  are  the  eigenvadues  cind  the  the  eigenvectors  of  the  matrix  [^4].. 

The  constants  Q  must  be  determined  from  the  boundary  conditions.  Note  that 
the  eigenvalues  and  eigenvectors  depend  on  the  frequency  u).  The  frequency  used 
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to  find  the  eigenvalues  and  eigenvectors  is  not  a  natural  frequency  of  the  system 
uidess  it  makes  the  determinant  of  the  boundary  conditions  vanish.  We  form  the 
boundary  condition  equations  by  applying  the  boundary  conditions  to  the  general 
solution  (6.37). 

As  an  example,  we  take  a  fixed-free  rod  for  which  the  boundary  conditions  are 
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The  resulting  matrix  equation  for  the  boundary  conditions  is 


(6.38) 
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^51 

(l>62e^^  <i>63e^^^ 

where  the  (f)ij  are  elements  of  the  eigenvectors  (the  second  index  j  indicates  the 
eigenvector  and  the  first  index  i  indicates  the  row  of  that  element),  and 
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Ni  =  {2as(f>i2  +  ‘2-o.Q^iz  +  4:OC3<l>i4)  (no  sum  on  z). 


The  natural  frequencies  of  extensional  vibration  axe  those  values  of  u  that  make 
the  determinant  of  the  matrix  in  (6.39)  vanish  (any  other  solution  to  the  above 
equation  is  the  trivial  one).  They  can  always  be  found  by  plotting  the  determinant 
as  a  function  of  w,  though  more  efficient  methods  can  be  devised.  Once  the  desired 
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natural  frequencies  are  found,  the  mode  shapes  are  determined  by  substituting  the 
corresponding  natural  frequency  back  into  (6.39),  cissuming  a  value  for  one  of  the 
coefficients  Q,  solving  (6.39)  for  the  remaining  coefficients  and  substituting  these 
coefficients  back  into  (6.37). 

In  order  to  create  some  plots  of  these  modes,  we  chose  a  steel  rod  and  fixed  the 
geometry.  Recalling  (2.52),  (5.13)  and  (5.14),  the  extensional  constitutive  coefficients 
and  inertial  coefficients  for  a  rectangular  rod  are 
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ufibh  libh^ 

,  O!io  =  <^11 
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12 


2u 


b^ 


12 


nhb^ 
12  ’ 


and  V  =  — ■ 
^  12 


The  equations  (6.39)  are  somewhat  ill-conditioned  for  structural  materials,  so  we 
use  a  relatively  thick  rod  with  L  =  0.1m,  b  =  0.02m  and  h  =  0.01m.  The  first 
four  natural  frequencies  for  this  rod  are  85.8,  256,  424  and  584  thousand  radians 
per  second,  and  the  mode  shapes  are  plotted  in  Figures  6.6  and  6.7.  The  first  three 
natural  frequencies  correspond  relatively  well  to  those  of  the  elementary  theory  for 
longitudinal  vibrations  of  a  bar  (the  theory  having  a  single  displacement  variable), 
but  later  modes  have  significant  interaction  with  the  lateral  extension.  Thus,  at  least 
for  a  thick  rod,  the  elementary  theory  diverges  from  a  more  complete  theory  after  the 
first  few  modes. 


64 


0  0.2  0.4  0.6  0.8  1 

4/L 

Figure  6.6:  The  first  and  second  extensional  modes  of  vibration  for  a  fixed-free 
(cantilevered)  Cosserat  curve  with  /x  =  82.7GPa,  =  7500Kg/m^,  u  =  0.30,  L  = 

0.1m,  b  =  0.01m  and  h  =  0.02m. 
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Figure  6.7:  The  third  and  fourth  extensional  modes  of  vibration  for  a  fixed-free 
Cosserat  curve  with  /x  =  82.7GPa,  =  7500Kg/m^,  v  =  0.30,  L  =  0.1m,  h  —  0.02m 

and  h  =  0.01m. 
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Chapter  7 


Examples  in  Moderate  Rotation  Theory 

We  present  examples  here  to  show  how  the  theory  established  in  Chapter  5  can 
be  used.  In  order  to  insure  that  our  examples  are  valid,  we  will  have  to  choose 
appropriate  loads  and  lengths  so  that  the  rotation  and  strain  are  of  the  proper  order, 
as  we  are  now  dealing  with  a  nonhnear  theory.  It  will  become  apparent  that  the 
moderate  rotation  theory  is  sensitive  in  this  respect.  We  will  show  that,  while  a 
certain  geometry  is  acceptable  for  a  given  deformation,  it  is  not  valid  for  another 
deformation.  If  necessary,  the  theory  can  be  adjusted  by  changing  the  assumptions 
to  suit  a  desired  application,  but  we  do  not  pursue  this  issue  here. 

The  examples  are  again  restricted  to  the  case  of  initially  straight,  uniform  rods, 
and,  since  they  all  involve  cantilevered  rods,  we  again  use  the  modified  auxihary 
motion  with  S  =  I,  f  =0,  c  =  0  and  s(t)  =  0  in  (2.69).  We  shall  also  drop  the  tildes 
associated  with  the  modified  auxiliary  motion  in  this  chapter. 
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7.1  Static  Response  to  a  Uniform  Distributed 


Load 

As  a  first  example  of  moderate  rotation,  consider  the  case  of  a  distributed  load  in 
the  — direction  which  is  408  times  the  weight  of  the  rod.  Use  of  a  circulax  rod 
implies,  from  the  balance  laws,  that  the  displacements  ^12  =  ^21  =  ^i3  =  ^31  =  O’  This 
results  in  an  uncoupling  of  the  flexural  equations  and  leaves  them  linear.  Thus,  the 
flexural  equations  are  solved  in  the  usual  manner,  and  the  solution  of  the  extensional 
equation  is  then  dependent  on  these  results.  The  parameters  wiU  be  carefully  chosen 
so  as  to  satisfy  the  assumptions  of  the  theory.  We  first  set  the  measure  Cq  =  0.0001. 
We  choose  steel  as  the  material  with  fi  =  82.7GPa,  pj  =  7500Kg/m^  and  1/  =  0.30. 
Next,  we  arbitrarily  choose  the  length  of  the  rod  as  one  meter,  but  the  applied  load 
is  selected  based  on  the  knowledge  that  it  would  produce  the  maximum  deformation 
623  =  ^0  at  the  fixed  end.  Finally,  the  radius  is  chosen  as  i2  =  0.1m  to  maJce  the  tip 
rotation  moderate.  The  assumptions  of  the  theory  regarding  the  partial  derivatives 
of  strain  and  displacement  will  also  be  satisfied  in  this  example,  as  will  be  shown  in 
Figure  7.1. 

7.1.1  Flexural  Response 

We  begin  by  solving  the  flexural  equations,  which  are 

(^23  +  ^32)  +  A/^  =  0,  (7.1) 

OlS  —  0^5  (^23  +  ^32)  =  0.  C^’^) 
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The  boundary  conditions  for  flexure  axe,  again, 


U2(0)  =  0,  ^23(0)  =  0,  —  as  (S^siL)  +  ^32(-t'))  —  0, 

m^^{L)  =  ai5^^{L)  =  0.  (7.3) 

The  solution  is  obtained  by  first  integrating  (7.1)  with  respect  to  (  and  then  applying 
the  third  boundary  condition: 


05  (^23  +  ^32)  +  a  ~  L)  —  0. 


(7.4) 


Next,  we  substitute  this  result  into  (7.2),  integrate  and  apply  the  fourth  boundary 
condition: 

r.2\ 

(7.6) 
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Integrating  again  and  applying  the  second  boundary  condition  gives 


^23  = 


Ar  te 


016  V  6 


2  2 


) 


(7.6) 


We  can  now  substitute  this  result  back  into  (7.4),  and  then  introduce  the  second 
of  (5.25)  to  recover  the  displacement  U2: 


as  as  ais  \  6  2  2  ^ 

Finally,  we  integrate  once  more  and  apply  the  first  boundary  condition: 


(7.7) 


U2  =  -A/^ 


as  V  2  j 


1  (e  Le  L^e' 

ais  V24  6  4 


(7.8) 
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7.1.2  Extensional  Response 


We  now  address  the  extensional  deformation.  Recalling  (5.20),  the  equihbrium  equa¬ 


tion  is 


—  ^4q;3  ^^33  4-  -^32^32^  + 


'^s^ 


23 


The  boundary  conditions  for  extension  are 


di 


=  0. 


(7.9) 


'a^23. 


„3(0)=0,  n^(L)=4a3(M^)  +  ^M^)M^))+«i5(^^Wj  =  »•  (7-10) 

The  second  boundary  condition  reveals  that  £*33  =  0  throughout  the  rod.  We  integrate 
the  equilibrium  equation,  recall  the  third  of  (5.25)  and  apply  the  second  boundary 
condition  to  get  the  simple  result 


-  5u3  1^  ^  ai5 

033  =  ^7“  =  --O32O32  7— 

oi  2  4^3 


(t)'- 


(7.11) 


The  result  U3  is  obtained  by  integrating  (7.11),  and  the  constant  of  integration  van¬ 
ishes  due  to  boundary  condition  (7.10)i: 


U3 


Next,  we  plot  the  displacements.  First,  it  should  be  noted  that  the  coefficients  03, 
0:5  and  ais  are  given  in  Section  5.2.  The  applied  load  is  Af*  =  —408gpQT)2  = 
—4000 Pq  D2,  which  yields  Xp  =  —4000 pQ  ir  B?  and  =  0.  The  numerical  results  are 
shown  in  Figures  7.1  and  7.2,  where  we  have  also  plotted  the  results  using  the  initial 
strain-displacement  relations  from  the  alternate  definition  of  moderate  rotation  (4.17) 
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(corresponding  to  Casey  and  Naghdi  [6,  §3.3])  to  show  that  the  two  methods  are  in 
agreement  to  the  order  of  approximation  prior  to  making  any  simphfication. 

Reviewing  the  component  results  in  Figure  7.1,  it  is  clear  that  the  supremums  of 
the  components  S23  and  S32  increase  only  by  a  factor  of  three  upon  taking  their  partial 
derivatives,  while  the  components  themselves  are  clearly  moderate.  The  measure  623 
and  its  partial  derivative  with  respect  to  ^  are  also  seen  to  be  small,  as  required  by 
the  assumptions  of  the  theory.  We  therefore  deem  the  theory  to  be  valid  for  this 
example. 

7.2  Free  Flexural  Vibration  in  the  First  Mode 

We  consider  here  the  fundamental  flexural  vibration  of  the  rod  used  in  Example  1.  The 
beam  will  be  released  from  rest  in  the  mode  one  position  (it  is  understood  that  both 
U2(^)  and  A23(0  must  begin  in  the  fundamental  mode),  and  will  therefore  continue 
in  mode-one  flexure  perpetually.  In  the  moderate  rotation  theory  this  will  cause  a 
forcing  input  to  the  extensional  vibration  equations  that  wiU  result  in  an  extensional 
vibration.  To  make  the  theory  valid,  we  need  only  verify  that  the  ampHtude  satisfies 
the  assumptions  of  the  moderate  rotation  theory. 

7.2.1  Flexural  Response 

We  recall  first  that,  in  the  moderate  rotation  theory,  the  flexural  equations  are  in¬ 
dependent  of  the  extension  and  are,  in  fact,  the  linear  flexural  equations.  Thus, 
we  may  apply  the  results  of  the  infinitesimal  theory  but  with  a  vibration  satisfying 
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Figure  7.1:  Static  deformations  of  a  cantilever  rod  due  to  load  Xp 
with  fx  =  82.7GPa,  =  7500Kg/m^,  i/  —  .3,  L  =  1.0m  and  R  =  0.1m. 
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Figure  7.2:  Static  displacements  of  a  cantilever  rod  due  to  load  A/^  =  — 4000/?o  A 
with  /X  =  82.7GPa,  Po  =  7500Kg/in^,  u  =  .3,  L  =  1.0m  and  R  =  0.1m. 


the  assumptions  of  the  moderate  rotation  theory.  Using  the  geometry  and  mate¬ 
rial  of  the  example  of  section  6.1,  we  iterate  the  frequency  equation  to  find  that 
(jip  =  919.59rad/sec.  We  substitute  this  frequency  into  the  mode  shape  and  seek  an 
amplitude  Ai  that  makes  the  root  strain  small  but  the  tip  rotation  moderate.  Letting 
€q  =  0.001,  an  amplitude  of  Ai  =  1/70  approximately  satisfies  these  requirements. 
The  deformation  U2{0  is  shown  in  Figure  7.3,  which  also  shows  several  of  the  other 
deformation  measures  needed  for  future  analysis.  One  can  see  clearly  from  this  fig¬ 
ure  that  the  supremum  of  2e23  is  of  O(eo)  while  the  supremum  of  A32  is  moderate 
(0(ey^)).  What  is  interesting  about  this  figure  is  that  we  can  use  it  to  observe  the 
order  of  all  terms  in  the  balance  of  director  momentum  equation.  Upon  doing  this  we 
discovered  that  while  (after  dividing  through  by  as)  both  terms  on  the  left-hand-side 
of  (6.2)  are  of  O(eo),  the  right-hand-side  (the  difference  between  the  other  terms)  is 
smaller  than  This  circumstance  is  not  expected  from  the  theory  nor  is  it 

excluded  by  the  assumptions. 


7.2.2  Extensional  Response 


We  now  consider  the  extensional  vibration.  This  is  governed  by 


dSss  d^us  8632  7  r,  8S32 

4a3-;^ - =  -403-^032  “  2ai5- 


dt^ 


’ 


(7.13) 


where  ^32  =  ^sin(a;t  -  (/>)  and  S23  —  A23sin(a;t  —  </>).  Since  the  right-hand  side 
of  this  equation  is  a  known  function  of  ^  and  t,  it  can  be  solved  as  an  uncoupled 
partial  differential  equation.  Strictly  speaking,  the  initial  conditions  of  the  extensional 
deformation  must  also  be  considered;  however,  for  our  purposes  such  consideration 
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Figure  7.3:  The  flexural  deformations  due  to  the  first  mode  of  a  cantilevered  rod 
having  small  strain  and  moderate  rotation  with  =  82.7GPa,  Pq  =  7500Kg/m^, 


=  0.30.  L  =  1.0m  and  R  =  0.1m. 


would  take  us  far  beyond  the  point  of  this  example  and  would  only  involve  standard 
methods  (cf.  [42,  Ch.  7]).  Therefore,  we  shall  consider  the  initial  conditions  of  the 
extensional  deformation  to  be  such  that  only  the  particulax  solution  (also  known  as  the 
steady-state  response  in  the  vibration  literature)  survives.  Without  loss  in  generality, 
we  set  (f)  =  0  (i.e.,  S32  =  dU^ld^  sinuit  and  S23  =  A23sinct;t,  an  assumption  which 
only  involves  a  shift  in  time  of  the  harmonic  response)  and  rec2dl  the  second  and  third 
of  (5.25)  to  transform  the  above  equation: 


^  ^  d^U3 


f  d?U2dU2  d^A23dA23\  .2  ^  fn  i  a\ 


„  i‘Ui  dUt 


d?'  A23  dA23 


d^2  di 


(1  -  cos  (2wij,t)) . 


The  appropriate  boundary  conditions  axe 

W3(0,t)  =  0,  n^(L,t)  =  4a3  2  =  0. 

(7.15) 

The  forcing  frequency  for  extensional  vibration  is  tAvice  the  vibration  frequency  of  the 
flexural  response.  There  is  also  a  static  component  to  the  forcing  function,  which  is 
clearly  apparent  after  noting  the  right-hand-side  of  (7.15). 

The  solution  is  compHcated  by  the  second  boundary  condition,  which  is  time  de¬ 
pendent.  Meirovitch  [42,  pp.  300-308]  discusses  a  method  of  solving  such  equations. 
We  win  take  his  approach  but  provide  a  minimum  of  details  on  the  method.  The  ap¬ 
proach,  which  he  partially  attributes  to  a  statement  made  by  Courant  and  Hilbert  [9, 
p.  277],  involves  assuming  a  solution  of  the  form 


=  +  (7.16) 
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where  the  function 


f{t)  =  4a3^{L,t)  =  -2a3  “  "is  > 


(7.17) 


7  \d( 

is  introduced  to  account  for  the  time-dependent  boundary  condition,  and  v{(,t)  is 
a  function  intended  to  have  only  homogeneous  boundary  conditions.  The  function 
h{^)  must  then  be  established  to  satisfy  the  boundary  conditions.  In  terms  of  these 
functions,  the  balance  law  and  boundary  conditions  are 


d^v 


d'^v . 


=  -  (  2a3 


(PUo  dUo 


de  d( 


+  q;i5- 


dP'  A23  d^23 


de  di 


(1  —  cos  {2uJipt))  — 


<ie 


iP  ' 


(7.18) 


n(0,t)  =  -h(0)/(t),  ia3^{L,t)  =  /(<)  ^1  -  4a3^{i:))  •  (7.19) 

To  render  the  boundary  conditions  on  the  function  v{^,t)  homogeneous,  the  condi¬ 
tions  on  the  function  h{()  must  be 


h{0)  =  0,  4a3^{L)  =  1. 


(7.20) 


There  are  many  simple  functions  that  satisfy  these  conditions.  In  a  similar  example, 
Meirovitch  [42,  §7-14]  used  a  step  function.  Here  we  chose  the  smooth  function 

L 


'“<«  =  2^5:  (1  - (^f)) 


(7.21) 


We  are  now  in  a  position  to  find  the  response 

We  begin  by  finding  the  natural  frequencies  and  normalized  eigenfunctions.  This 
is  done  in  the  usual  way,  separating  variables  and  applying  the  boundary  conditions 
to  find 
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We  now  introduce  the  modal  coordinates  T]n{t)  by  invoking  the  expansion  theorem 
(i.e.,  assuming  the  response  t)  to  any  excitation  can  be  represented  by  an  infinite 
series  in  the  eigenfunctions): 


Next,  we  substitute  this  assumed  solution  into  (7.18)  to  get 


(2n  —  1)7 

2r 


2L 

'  / 

/r>  d^U2  dU2  (P A23  dA.23 

-  r’le  IT 

>2  ■1‘m 


(7.23) 


j  (1  -  cos  (2ui,t))  - 


(7.24) 


We  now  take  advantage  of  the  orthogonality  of  the  eigenvectors  V„(()  by  mnltiplying 
this  entire  equation  by  Kn(^)  and  integrating  along  the  length  of  the  rod.  In  doing  so, 
only  the  m*'*  terms  of  the  series  survive,  and  the  resulting  modal  equations  of  motion 
become  uncoupled: 

2  /J.\  Tvr  /  \ 

—^^2  +^mnm{t)  =  (7.25) 

where 
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where 


Vm  =  —  ,  ^  (cos  7r(m  -  1)  -  1)  4-  —  (cos  mTT  -  1)  ,7^  =  1/^+  .  .• 

27r  [7r(m  - 1)  ^  ^  ^  '  m  '  'J  7r(2m  - 1) 


Note  that  the  first  term  in  the  expression  for  Vm  is  understood  to  vanish  for  m  =  1. 

The  responses  r}m{i)  can  be  resolved  into  the  static  and  dynamic  components 
respectively  as  7/m(f)  =  'Hms  +  One  can  immediately  solve  for  the  static 

components: 


(2m  -  1)2 


-  [^\f 
7r2V  2  [Jo 


AT  [  .  /(2m-J> 


(PU2dU2  ai5d2A23dA23\  ,, 

dC  ^2a3  de  de  j  ^ 


■KUm  ( dUi,  . 

4L  [  di 


(7.27) 


For  the  dynamic  part,  we  assume  a  solution  of  the  form  —  C'm  cos  2a;ipt, 

substitute  into  (7.25)  and  solve  for  the  coefficients  Cm- 


1  |■2n;2  ( fdU2,,y  .  a,,  fdA23,ry\  2 


a?2  —  4u)? 


(7.28) 


The  total  modal  responses  are  then  given  by 


“  ,.,2  _  4,„2 


1  2a;2 

■  4u>?  TT  i  I  d£ 


^2  —  4cc>i 


(7.29) 


Finally,  the  extensional  response  is  determined  by  substituting  these  modal  responses 
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back  into  the  series  (7.23)  and  substituting  the  resulting  expression  into  (7.16)  to  get 


OO 

=  E  VMM)  - 

m—X 


where  we  have  separated  the  response  into  time  and  spatial  dependencies:  U3{(,t)  = 
U3^{^,t)  +  U3s{C)  =  U3j{^)cosuJipt  +  U3g{().  Note  that  there  is  a  “static”  displacement 
over  which  the  time  dependent  part  gets  superimposed.  Also  note  that  because  the 
natural  frequencies  of  the  extensional  modes  are  much  higher  than  the  frequency  of 
the  flexural  vibration  (the  first  mode  is  a  factor  of  five  larger),  the  difference  between 
these  two  responses  is  quite  small.  This  makes  the  total  response  nearly  zero  at  the 
same  instant  of  time  that  the  flexural  response  is  zero.  We  note  that  our  choice  for 
h{^)  results  in  rapid  convergence  of  the  series,  and  we  show  the  sums  of  the  first  eight 
terms  in  Figure  7.4.^ 

^[42]  contains  a  footnote  stating  that  many  different  forms  for  h{^)  may  be  suitable  in  a  given 
problem.  We  know  of  no  proof  that  all  such  functions  converge  to  the  same  result,  but  in  using  the 
step  function  that  Meirovitch  used  in  his  example  [42,  Example  7.2]  we  achieved  the  same  result  as 
with  our  choice  of  h{^)  (7.21)  except  that  about  70  terms  were  required  to  reach  the  same  level  of 
convergence. 
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Figure  7.4:  Extensional  deformation  due  to  the  first  mode  of  flexural  vibration  in  a 
rod  with  jx  =  82.7GPa,  pQ  =  7500Kg/m^,  v  =  0.30,  L  =  1.0m  and  R  =  0.1m. 


7.3  PVee  Flexural  Vibration  in  the  Fifth  Mode 

We  consider  here  the  fifth  mode  of  flexural  vibration  of  the  same  rod  as  the  previous 
examples,  but  this  time  in  the  fifth  mode.  We  -will  show  that  the  assumptions  of  the 
moderate  rotation  theory  are  not  satisfied  for  this  motion.  To  make  this  clear,  we 
state  the  primary  assumptions  in  clear  terms  (using  our  alternate  definition  of  small 
strain  with  moderate  rotation)  as 

623  =  2  (^23  +  ^32)  =  ^(^0)  Co  0, 

1023  =  -  (A23  —  A32)  =  O(eo )  as  60  0,  (7.32) 
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where  we  are  now  considering  the  A.y  to  be  the  values  of  the  corresponding  Sij  at  the 
time  of  majcimum  excursion  of  the  vibration.  We  seek  an  amplitude  As  for  the  fifth 
eigenfunction  (6.19)  that  satisfies  the  above  conditions  for  say  co  =  0.001  (cq'  «  0.03). 
We  show  the  numerical  values  of  the  deformations  A23  s-nid  ^32  for  an  amplitude  As 
that  satisfies  (7.31)  in  Figure  7.5,  where  we  also  show  the  numerical  results  for  the 
maximum  values  of  W23  and  623.  The  strain  and  rotation  at  this  amplitude  are  both 
small.  If  we  increase  the  amplitude,  both  measures  will  increase  by  the  same  factor. 
Thus  the  strain  and  rotation  are  always  comparable  in  order  (when  the  strain  is 
small,  the  rotation  is  also  small;  and  when  the  rotation  is  moderate,  the  strain  is  also 
moderate).  Consequently,  the  moderate  rotation  theory  cannot  be  applied  in  this 
case.  If  we  desire  to  study  the  theory  in  mode-five  flexural  vibration,  we  must  alter 
the  rod  to  suit  the  assumptions  of  the  theory.  Reduction  of  the  rod  radius  tends  to 
make  the  order  of  the  rotation  larger  than  that  of  the  strain.  If  we  reduce  the  radius 
of  the  rod  from  0.1m  to  0.02m,  the  mode-five  deformations  shown  in  Figure  7.6  result. 
In  this  case,  the  strain  can  be  small  while  the  rotation  is  moderate,  and  the  theory  is 
valid. 
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Figure  7.5:  Tlie  fifth  mode  of  flexural  vibration  for  an  infinitesimal  theory  with 
H  =  82.7GPa,  pl  =  7500Kg/m^  v  =  0.30,  L  =  1.0m  and  R  =  0.1m. 


0.04 


Figure  7.6:  The  fifth  mode  of  flexural  vibration  for  a  small  strain  and  moderate 
rotation  theory  with  //  =  82.7GPa,  =  7500Kg/m^,  v  =  0.30,  L  =  1.0m  and 

R  =  0.1m. 
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Chapter  8 


Free  Vibration  of  a  Whirling  Rod 


8.1  Introduction 


The  study  of  the  free  vibrations  of  a  rod  whirhng  about  a  fixed  end  has  a  long  history. 
Here,  the  rod  is  spinning  about  a  fixed  axis,  as  indicated  in  Figure  8.1.*  The  whirling 


Figure  8.1:  A  Cosserat  curve  whirling  about  the  Di(0)  axis. 


of  the  rod  alone  induces  stretching  in  an  unconstrained  theory  and  may  induce  other 

*This  figure  also  illustrates  various  bases,  which  we  will  subsequently  define.  Also  note  that 
S  =  S(t)  is  conveniently  taken  to  be  the  rotation  tensor  of  the  shaft  to  which  the  rod  is  attached. 
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deformations  as  well,  depending  on  both  the  material  and  geometric  symmetries  of 
the  rod.  In  addition  to  these  deformations,  the  rotation  also  affects  any  superposed 
vibrations,  altering  both  the  frequencies  and  shapes  of  the  modes  as  compared  to  a 
stationary  cantilevered  rod.  Due  to  the  variety  of  applications,  it  is  of  considerable 
interest  to  examine  these  vibrations  to  determine  the  spectral  responses  and  linear 
stabihty  of  the  whirling  rod. 

This  problem  interests  us  not  only  because  of  its  apphcation  to  turbomachinery, 
propellers  and  hehcopter  rotors,  but  also  because  our  properly  invariant  auxiliary 
motion  is  well  suited  to  this  type  of  problem.  While  aerodynamic  forces  are  generally 
important  in  applications  and  can  easily  be  accounted  for  in  the  Cosserat  theory  (in 
the  terms  Af  and  Al“),  we  do  not  consider  them  in  this  chapter.  This  chapter  serves 
as  an  example  of  how  O’Reilly’s  [50]  properly  invariant  approximate  theory  can  be 
used  and  provides  a  contrast  to  the  usual  approach  talcen  in  approximate  rod  theories. 


8.2  Previous  Work  on  Whirling  Rods 

The  hterature  abounds  with  work  on  whirhng  rods.  Leissa,  in  his  review  article  [38], 
provides  102  references  pertaining  only  to  rotating  blades,  and  many  more  articles 
have  appeared  since  his  review.  Another  review  article,  by  Rao  [53],  lists  140  refer¬ 
ences.  Most  of  these  articles  address  the  vibratory  response  under  constant  rotation 
speed,  but  numerous  effects  have  been  repeatedly  considered.  These  include  flexu¬ 
ral,  extensional^  and  torsional  vibrations;  both  coupled  and  uncoupled.  They  also 

^By  extension,  we  mean  here  longitudinal  extension  (see,  for  example,  Anderson  [1]  or  Hodges  and 
Bless  [29]).  The  lateral  deformations  (due  to  the  Poisson  effect)  are  not  discussed  in  the  references 
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include  the  effects  of  hub  radius,  attachment  angle,  tapering,  pretwist,  asymmetry 
of  the  reference  configuration,  and  material  nonlinearity.  However,  the  rod  models 
used  are  generally  established  in  an  od  hoc  manner  using  variational  principles  com¬ 
bined  with  kinematical  and  constitutive  assumptions  (see,  for  example,  [1],  [3]  and 
[37]).  Furthermore,  the  correct  formulation  of  the  problem  involves  delicate  scaling 
arguments.  Surprisingly,  there  is  still  controversy  in  the  literature  on  these  issues 
(cf.  the  discussions  of  Hodges  [30]  and  Ko  [36],  for  example).  There  is  also  a  large 
body  of  related  research  in  the  astronautical  literature  (see,  for  example,  Levinson 
and  Kane  [35],  [39]  and  Simo  and  Vu-Quoc  [58],  [59]). 

An  improved  rod  model  will  help  to  solidify  the  theory  of  whirhng  rods.  For 
example,  Bhuta  and  Jones  [2]  used  the  one- dimensional  wave  equation  to  study  the 
effects  of  rod  rotation  on  extensional  vibrations  and  predicted  a  decrease  in  the  nat¬ 
ural  frequencies  due  to  the  rotation.  Hodges  [28]  later  pointed  out  that  the  natural 
frequencies  of  the  extensional  vibrations  may  actually  increase  rather  than  decrease 
compared  to  the  fixed  rod,  depending  on  the  material  nonlinearities,  however,  his 
model  was  also  highly  simplified  (we  have  seen,  in  Chapter  6,  that  lateral  extensions 
can  be  included  by  use  of  the  Cosserat  curve  model).  His  discussion  should,  however, 
caution  us  about  drawing  general  conclusions  regarding  the  effects  of  the  rotation  on 
the  eigenstructure  of  a  given  rod.^ 

Another  issue  of  controversy  concerns  the  inertial  terms.  Most  previous  researchers 

(for  example,  [l],  [2],  [28],  [62]  and  [37])  have  chosen  to  neglect  the  Coriolis  terms 
we  cite  in  this  section, 

^See  Hodges  [27],  Venkatesan  and  Nagaraj  [62]  and  Hodges  and  Bless  [29]  for  further  details  on 
this  effect. 
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coupling  the  exteusional  and  flexural  vibrations  of  an  otherwise  uncoupled  straight 
rod.  After  neglecting  the  Coriolis  terms,  Lee  and  Lin  [37]  recently  obtained  exact 
flexural  eigensolutions  for  the  cases  where  the  coefficients  in  the  resulting  partial 
differential  equations  are  polynomials.  Rcio  and  Carnegie  [54]  and  [55]  did  consider 
Coriolis  effects,  but  they  used  the  nonhnear  theory  derived  by  Carnegie  [3],  which 
differs  substantially  from  other  theories,  including  ours.^  It  is  our  opinion  that  the 
Coriolis  terms  should  not  be  neglected  in  an  eigensolution  unless  it  is  restricted  to 
cases  where  applied  forces  sufficiently  dominate  such  terms.  It  is  interesting  to  note 
that  Leissa  [38]  points  out  that  there  is  stiU  debate  on  this  point. 

8.3  Properly  Invariant  Small  Deformations 
Superposed  on  a  Large  Deformation  of  a 
Whirling  Rod 

O’Reilly  [50],  while  not  specifically  addressing  the  whirling  rod,  provided  a  properly 
invariant  approximate  theory  of  rods  for  small  motions  superposed  on  a  large  motion 
that  is  well  suited  to  this  problem.  His  theory,  which  is  based  on  the  Cosserat  theory 
of  Green  and  Naghdi,  uses  momentum  balance  laws  to  determine  the  equations  of 
motion.  It  also  accounts  for  arbitrary  large  motions  accompanied  by  infinitesimal 
superposed  deformations.  The  resulting  equations  take  on  a  more  familiar  form  after 

being  written  in  terms  of  the  auxiliary  motion,  which  is  equivalent  to  referring  the 

^Carnegie  [3]  neglects  longitudincd  extension  but  considers  axial  motion  due  to  the  integrated 
effects  of  flexure,  which  results  in  the  presence  of  nonlinear  Coriolis  terms  in  his  equations. 
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motion  to  a  corotating  frame  except  that  it  is  properly  invariant  under  superposed 
rigid  body  motions. 

We  will  use  the  unconstrained  theory  in  our  derivation,  formulating  the  problem 
as  one  of  small  vibrations  superposed  on  a  large  steady  deformation.  The  large 
deformation  is  caused  by  the  rotation  of  the  rod  at  a  constant  angular  velocity.  Many 
of  the  subsequent  developments  parallel  those  of  Green,  Knops  and  Laws  [16],  and  we 
refer  the  reader  to  their  work  for  details,  especially  concerning  the  constitutive  laws 
for  a  theory  of  small  deformations  superposed  on  a  large  deformation.  For  simplicity, 
we  perform  most  of  the  analysis  on  a  corotating  reference  frame,  but,  in  the  case 
of  an  unconstrained  Cosserat  curve,  such  analysis  is  formally  equivalent  to  the  use 
of  the  properly  invariant  modified  auxiliary  motion.  Consequently,  our  results  will 
be  properly  invariant  under  superposed  rigid  body  motions  according  to  the  theory 
developed  by  O’Reilly  [50].  Because  we  do  not  anticipate  solving  for  the  large  steady 
deformation,  we  leave  the  rod  material  and  geometry  general.  In  particular,  we  permit 
the  possibility  that  the  rod  is  composed  of  an  anisotropic  material.  As  discussed  by 
Green  and  Naghdi  [23],  the  Cosserat  theory  they  developed  is  sufficiently  general  to 
encompass  this  case. 

By  a  steady,  large  deformation,  we  refer  to  the  deformation  which  exists  in  a 
rod  whirling  at  a  constant  angular  velocity.  As  we  cire  restricting  our  theory  to 
elastic  rods,  we  presume  that,  but  do  not  consider  the  mechanisms  by  which,  the 
transient  vibrations  (those  that  normally  accompany  a  spin-up  from  a  stationary  rod 
to  one  spinning  at  a  constcint  angulcir  velocity)  subside.  This  large  deformation  is 
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one  in  which  the  components  of  all  vector-valued  quantities,  when  referred  to  the 
corotational  basis  {SD,(0)},  axe  independent  of  time.  This  implies  that  the  properly 
invariant  strain  measures  associated  with  this  motion  are  independent  of  time.  We 
note  that  this  definition  is  consistent  with  the  notion  of  a  steady  motion  used  in 
the  literature  on  whirling  rods.  Once  the  large  motion  is  determined,  superposed 
vibrations  may  be  considered. 

The  superposed  small  deformations  we  consider  axe  infinitesimal  vibxations.  To 
obtain  the  balance  laws  for  these  deformations,  it  is  appropriate  to  hneaxize  the 
balance  laws  for  the  Cosserat  curve  about  the  large  deformation.  For  consistency 
with  the  notation  of  O’ReiUy  [50],  we  use  the  notation  2(  )  for  the  full  motion  and 
i(  )  for  the  large  motion.  We  also  use  the  designation  21  (  )  for  the  superposed  or 
difference  motion.  Thus,  the  vaxiables  of  the  motion  may  be  decomposed  as 

2n(^,t)  =  in(^,t)  -I-  2in(^,t),  2m“(^,t)  =  -f-  2im"(^,t),  (8.1) 

2k“(^,  t)  =  ik“(^,  t)  +  2ik“(e,  t),  2U(^,  t)  =  iu(^ ,  t)  +  2lU(^,  t),  (8.2) 

2^/3(^,<)=  2d/3(,f,t)  =  D;3(^)-t- 2^^(^,<)  =  ld^(^,  <)+ 2l^/3(^,  <). 

(8.3) 

We  note  here  that  the  boundary  conditions  for  the  problem  axe 

2U(0,  t)  =  0,  2da(0,  t)  =  S(t)Dc(0),  2n(i,  t)  =  2m“(T,  t)  =  0.  (8.4) 
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8.4  Balance  Laws  and  Responses 

We  begin  with  the  balance  laws  for  the  modified  auxihary  motion  associated  with 
the  full  motion  of  the  Cosserat  curve,  fixing  the  origin  of  our  reference  frame  at  the 
point  where  the  directed  material  curve  is  attached  to  the  rotating  shaft.  The  right- 
hand-sides  of  the  balance  laws  are  easily  obtained  from  the  counterparts  of  (2.67) 
and  (2.68)  for  the  modified  auxiliary  motion.  For  the  problem  of  interest,  we  choose 
S  =  S(t),  =  0,  c  =  0  and  s  =  0  in  (2.69).  The  non-trivial  balance  laws  are 

^  =  AS^(t)  [«''S(t)  2? (e,  t)  -t-  2fiS{t)  2v(6 t)  -f-  S(i)  2’i^(e,  t)+ 

/  2d^(^,  t)  +  2ns{t)  2W^(e,  t)  +  S{t)  <))]  ,  (8.5) 

-  2k“  =  AS^(t)  [2/"''(n2S(t)2d^(e,<)  +2fiSit)2WfiU,t)  +Sit)2-^fi{U))  + 
[n^S{t)  2r(^,  t)  +  2ns{t)  2V(^,  t)  -b  S{t)  2HC,  <))]  ,  (8-6) 

where  O  =  SS^  represents  the  angular  velocity  tensor,  which  is  assumed  to  be 
constant.  Note  that  the  form  of  the  balance  laws  are  such  that  y"  0. 

The  tilde  notation  in  (8.5)  and  (8.6)  refers  to  vectors  of  the  modified  auxiliary 
motion  (for  example,  dj}  =  S^d/3).  As  the  pivot  rotation  tensor  R(|^  =  0)  differs 
from  the  shaft  rotation  tensor  when  shear  is  present  at  the  cantilevered  end  of  the 
rod  (and  shear  generally  is  present,  even  for  the  steady  motion,  when  the  rod  lacks 
adequate  symmetry),  we  have  not  used  the  standard  auxihary  motion  in  this  problem. 
Furthermore,  the  rotation  tensor  at  the  pivot  due  to  shear  in  the  steady,  large  motion 
may  itself  be  large. 
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Assuming  the  referential  basis  of  the  rod  D,(^)  to  be  orthonormal  and  the  rotation 
of  the  shaft  to  be  about  the  Di(0)  axis,  the  tensors  S(t)  and  fi  are 

S{t)  =  Di(0)  (Si  Di(0)  +  cos  fit  (D2(0)  ®  D2(0)  +  D3(0)  (S)  03(0))  + 

sin  nt  (D2(0)  (Si  03(0)  -  D3(0)  (g)  D2(0)) ,  (8.7) 

fi  =  S(t)S^(t)  =  fi  (D2(0)  ®  D3(0)  -  D3(0)  (81  D2(0))  .  (8.8) 

The  following  identities  axe  also  noted: 

S^{t)f2^S{t)  =  -fi2  (D2(0)  ®  D2(0)  +  D3(0)  (Si  03(0)) ,  (8.9) 

2S^{t)nS{t)  =  =  2fi  (D2(0)  (SI  D3(0)  -  D3(0)  (8  D2(0))  .  (8.10) 

Now  note  that  2f(^>^)  =  1^(0  +  2u(^,t)  and  2d^(^)0  =  +  2^/3(^5i)>  so 

that  (8.5)  and  (8.6)  can  be  written  as 

^  =  AS^(t)  [«^S(t)(R(0+  2u(^,t))^-2^^S(<)2v(e,t)  +  S(^)2<^(^,t)+ 

/  (f72S(<)  (D^(0  +  2hiU))  +  2f2S{t)2WfliU)  +  S{t)2^0iU))] ,  (8.11) 

-  2k"  =  AS^(t)  [t/“^  [f2^S{t)  (pfiiO  +  2h{U))  +  2f2S{t)2^i3{U)+ 

S(t)2W/3(^,i))  +  i/“  (!7^S(t)  (R(0  +  2u(<^,<))  +217S(t)2v(^,t)+ 

S(t)2<^(^,t))].  (8.12) 

Considering  (8.11)  and  (8.12)  in  conjunction  with  (8.9),  it  is  clear  that  there  are 
time-independent  terms  on  the  right-hand-sides  of  (8.11)  and  (8.12)  which  produce 
the  large  deformation  due  to  the  steady  rotation.  We  shall  first  consider  this  large 
deformation  and  then  consider  superposed  vibrations.  As  we  have  written  (8.11)  and 
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(8.12)  in  terms  of  modified  auxili«iry  motion  variables,  it  should  be  clear  that  all  field 
variables  associated  with  the  first  modified  auxihary  motion  |  if,  id^  j-  are  functions 
of  ^  only  rather  than  being  functions  of  both  ^  and  t,  as  they  are  for  the  motion 

{ir,  ida}. 

8.4.1  Response  to  a  Steady  Rotation 

If  we  incorporate  the  decompositions  (8.1)  -  (8.3)  into  (8.11)  and  (8.12),  we  obtain 
the  balance  laws  associated  with  the  steady  rotation; 

[R(0  +  iu(0  +  /  (D/3(0  +  i^'/3(0)] 

=  AS^(<)(iV  +  /iW;i),  (8.13) 

ik“(0  =  AS^(t)«"S(t)  [3/“^(D^(0+  i^/3(0)  +  2/“(R(0+  i«(0)] 

=  AS^(t)  (j/“^  +  y°‘  iv)  .  (8.14) 

The  constitutive  laws  for  this  large,  steady  motion  are  given  by  (2.34)  -  (2.36),  but, 
as  the  explicit  form  of  the  free  energy  V*  is  unspecified,  there  is  no  benefit  at  this  stage 
to  directly  substituting  the  constitutive  laws  into  (8.13)  and  (8.14). 

We  henceforth  assume  that  a  solution  ifi,  iih",  ik",  lU  and  to  the  steady 
rotation  problem  (8.13),  (8.14)  and  (8.4)  exists,  and  turn  to  considering  vibrations 
superposed  on  this  steady  motion. 
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8.4.2  Superposed  Vibrations 

In  the  superposed  motion  of  the  Cosserat  curve,  we  assume  that  the  vectors  2in, 
2iih“  and  21k"  are  of  O(eo)  as  €0  0,  where  Iq  now  refers  to  the  supremum  of  21E  = 

2E  —  lE  =  2E  —  lE.  We  also  assume  the  displacements  and  director  displacements 
21U  and  2i^a  are  of  O(co)  as  cq  0.®  To  obtain  the  balance  laws  for  the  superposed 
motion,  we  subtract  the  balance  laws  for  the  steady  motion  (8.13)  and  (8.14)  from 
those  of  the  full  motion  (8.11)  and  (8.12)  to  obtain 

=  AS^  21U  +  20s  2ifi  +  S  2ifi  +  / 

di  L 

=  AS^  [( 2V  -  iv)  +  ( 2W^  -  iW;9)]  =  h,  (8.15) 


52iiii" 


2ik“  =  AS^  [y"  21U  +  2ns  2ifi  +  S  2ifi)  + 

nSfi  +  2f2S  21^  +  S  21^/3)] 


=  AS' 


y"  ( 2V  -  iv)  +  y"'^  ( 2W/3  -  iw^)]  =  p“. 


(8.16) 


To  obtain  the  corresponding  component  balance  laws,  we  resolve  the  kinetic  vec¬ 
tors  2ih,  2im“  and  2ik“  onto  the  basis  |id,-  =  S^id.  j  and  take  the  inner  products 
of  (8.15)  and  (8.16)  with  id'^: 

[(2V-  iv)  +  y^(2W^  -  1W/3)]  •  id^  =  (8.17) 

. -I-  21^“’  lA,-.^  -  21^"’  =  as’’  fy"  ( 2V  -  iv)  +  y“'^  ( 2^0  -  iW/3)]  •  id"'  = 

(8.18) 

where  i\J  =  id^  •  didi/d^. 

®  Clearly,  these  developments  parallel  those  for  the  development  of  an  infinitesimal  theory. 
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As  the  right-hand-sides  of  (8.17)  and  (8.18)  are  functions  of  only  the  deformation 
and  the  angular  speed  (which  we  have  assumed  to  be  constant),  these  equations 
represent  nine  coupled,  linear  (after  neglecting  terms  from  the  right-hand-side  that 
are  of  0{l\)  as  Cq  0),  homogeneous  partial  differential  equations  with  variable 
coefficients  subject  to  the  following  boundary  conditions:® 

2iit’(0)  =  21^1^(0)  —  21^2^(0)  =  21^1^ (0)  =  21  ^2^(0)  21^1  (O)  =  21^2  (®)  “ 

(8.19) 

„n‘(£)  =  nk'‘(L)  =  =  2.m*’(i)  =  2,A'"{£)  -  2,m“’(i)  =  0. 

(8.20) 

The  constitutive  laws  for  the  superposed  vibrations,  which  depend  on  the  large 
motion,  can  be  obtained  from  those  given  by  Green,  Knops  and  Laws  [16,  Eqs. 
(3.49)-(3.52)]  by,  among  others,  neglecting  all  temperature  dependent  terms  in  their 
equations.’^  In  addition,  due  to  our  use  of  the  modified  auxiliary  motion,  our  results 
are  properly  invariant  under  superposed  rigid  body  motions.  Upon  substitution  of 
the  constitutive  laws  and  the  known  solution  of  the  large  motion,  equations  (8.17) 
and  (8.18)  become  functions  of  the  three  displacements  =  2ifi  •  id*  and  the  six 
director  displacements  21  =  21 

To  solve  the  eigenvalue  problem,  we  propose  to  use  Galerkin’s  method  by  approx¬ 
imating  the  response  as  a  series  of  n  comparison  functions®  for  each  displacement 

®As  the  large  motion  already  satisfies  the  boundary  conditions,  they  need  only  be  applied  to  the 
superposed  pait  of  the  total  motion. 

^ Green,  Knops  and  Laws  consider  a  homogeneous  deformation  of  the  rod  as  the  large  deformation 
in  their  small-on-large  theory.  However,  the  extension  of  their  discussion  on  constitutive  equations 
to  the  non-homogeneous  large  deformation  considered  here  is  trivial. 

® Comparison  functions  must  satisfy  the  boundary  conditions  (8.19)  and  (8.20).  Clearly,  because 
of  this  requirement,  the  constitutive  laws  of  the  rod  must  be  substituted  into  the  natural  boundary 
conditions  (8.20)  before  comparison  functions  can  be  found. 
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vaxiable; 


=  2.^;’ =  E  «(')»({.-)*«)•  (8-21) 

Upon  substitution  of  the  constitutive  laws  and,  subsequently,  of  these  assumed  solu¬ 
tions,  each  component  balance  law  taJces  the  form 

tci  t,  +  E  E  <7;’  [E  =  0,  (8.22) 

t=l  Li=l  J  a=lj;=l  u=i  '  ■' 

where  C,-  and  are  linear  differential  operators.  Upon  multiplying  each  component 

balance  law  by  n  comparison  functions  and  integrating  the  resulting  equations  &om 
^  =  0  to  ^  =  X,  a  linear,  homogeneous  system  of  9n  second  order  ordinary  differential 
equations  of  the  following  form  results: 


-{-  G<^  -b  K<^  =  0.  (8.23) 

If,  for  the  system  of  partial  differential  equations  (8.17)  and  (8.18),  a  common  set  of 
comparison  functions  can  be  found,  the  matrix  M  will  be  symmetric,  the  matrix  G 
will  be  skew-symmetric,  but  the  matrix  K  may  be  neither.  The  peculiar  nature  of 
K  is  due  to  the  structure  of  the  linearized  equations  (8.17)  and  (8.18).  This  can  be 
observed,  in  part,  from  the  extensional  equations  (6.35).  If  different  sets  of  comparison 
functions  are  used  for  the  various  displacement  variables  2iw*(^,i)  and  the 

aforementioned  symmetries  and  skew-symmetries  will  be  lost,  as  each  equation  of  the 
form  (8.22)  will  only  be  multiplied  by  n  comparison  functions,  whereas  it  may  contain 
all  9n  comparison  functions.^ 

^As  each  of  the  nine  equations  (8.17)  and  (8.18),  which  are  of  the  form  (8.22),  must  be  weighted 
by  only  n  comparison  functions  prior  to  integration  in  Galerkin’s  method,  the  only  way  to  ensure 
that  each  term  is  weighted  by  its  own  comparison  functions  is  to  have  a  common  set  of  n  compar¬ 
ison  functions  that  satisfy  all  of  the  boundary  conditions.  Finding  such  a  set  may  be  difficult,  as 
the  natural  boundary  conditions  (especi^Jly  ti^(L)  =  0)  may  then  involve  the  sum  of  compcirison 
functions  and  their  derivatives. 
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Due  to  the  lack  of  symmetry  of  K,  efficient  solution  techniques,  such  as  the  method 
developed  by  Meirovitch  [43]  and  used  by  Wickert  and  Mote  [65]  cannot  be  employed. 
Although,  in  some  cases,  one  may  be  able  to  exploit  the  symmetries  in  the  problem, 
one  must  usually  expect  to  use  a  general  complex  eigensolution  technique,  such  as  the 
complex  Lanczos  or  matrix  iteration  methods.  The  eigensolution  of  the  approximate 
system  yields  n  approximate  natural  frequencies.  The  approximate  eigenfunctions 
are  obtained  by  back  substitution  of  the  eigenvectors  into  (8.21). 

8.5  Discussion  and  Conclusions 

Although  the  approach  we  have  taken  is  much  easier  to  follow  than  many  of  the 
derivations  available  in  the  literature,  it  is  difficult  to  draw  general  conclusions  re¬ 
garding  the  natural  frequencies  of  the  whirling  rod  as  compared  to  those  of  the  fixed 
rod.  If  the  steady  motion  is  a  large  motion,  as  most  researchers  assume,  then  consti¬ 
tutive  laws  for  large  motions  must  be  used,  and  these  are  gener^dly  unavailable.  The 
difficulty  in  drawing  general  conclusions  is  compounded  by  the  complex  materials  and 
geometries  used  in  many  applications.  We  have  kept  our  theory  general  so  that  it  can 
be  used  for  all  of  the  important  applications. 

In  closing,  we  note  that  our  development  has  included  the  Coriolis  accelerations  in 
a  consistent  manner.  In  the  equations  for  the  vibrations,  these  accelerations  result  in 

a  non-trivial  G,  and  significantly  alter  the  eigenfrequencies  and  eigenfunctions  of  the 

^°The  only  generalization  we  have  omitted  is  that  of  allowing  for  a  hub  radius,  but  this  can  be 
easily  incorporated. 
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rod.^^  Finally,  our  formulation  is  unique  in  its  incorporation  of  lateral  deformations 
due  to  the  steady  motion. 


^^See  Wickert  and  Mote  [65]  for  examples. 
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Chapter  9 


Conclusions 

We  have  analyzed  four  approximate  rod  theories,  namely  the  infinitesimal  theory,  the 
theory  of  small  strain  accompanied  by  moderate  rotation,  the  theory  of  moderate 
strain  accompanied  by  small  rotation  and  the  theory  of  moderate  strain  and  rota¬ 
tion.  Two  different  approaches  were  given.  Also,  several  important  applications  of 
a  Cosserat  curve  model  were  analysed.  We  have  placed  particular  emphasis  on  the 
theory  in  which  small  strain  is  accompanied  by  moderate  rotation.  In  this  theory,  we 
constrained  the  lateral  extensions  to  illuminate  certain  features  of  the  balance  laws. 
In  general,  the  linear  balance  and  constitutive  laws  do  not  apply  when  the  strain  or 
rotation  is  moderate.  To  obviate  this  issue,  our  theory  and  illustrations  axe  based  on 
a  specific  quadratic  form  of  the  free  energy. 

In  the  development  of  the  constrained  theory  of  small  strain  accompanied  by  mod¬ 
erate  rotation,  we  improved  the  invariant  theory  of  O’Reilly  [50].  We  showed  that  the 
assigned  forces  and  assigned  director  forces  in  a  superposed  motion  need  not  be  ob¬ 
jective  unless  the  Lagrange  multipliers  are  objective.  We  argued  on  physical  grounds 
that  such  objectivity  would  be  too  restrictive,  and  referred  to  counter-examples  given 
in  O’ReiUy  and  Turcotte  [52]. 

The  first  two  applications  we  analyzed  pertained  to  the  flexural  and  extensional 
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vibrations  using  the  infinitesimal  theory.  As  the  hnear  equations  for  flexural  vibration 
of  a  Cosserat  curve  axe  equivalent  to  those  of  the  Timoshenko  beam  theory,  there  are 
substantial  results  available  in  the  Hterature,  but  we  proved  the  existence  of  previously 
undiscovered  modes.  In  solving  for  the  eigenfrequencies  and  eigenfunctions  of  the 
extensional  vibration,  we  have  found  the  first  such  results,  for  a  non-circulax  rod, 
that  include  lateral  extensions. 

We  gave  three  examples  in  the  moderate  rotation  theory.  In  the  first  of  these,  a 
static  response  to  a  distributed  load,  we  emphasized  the  fact  that  the  rod  geometry 
cind  loading  must  be  chosen  carefully  in  order  for  the  theory  to  be  valid.  The  second 
example,  that  of  flexural  vibration  in  the  first  mode,  emphasized  the  approach  that 
must  be  taken  to  determine  the  extensional  response  as  the  flexural  deformations 
serve  as  excitations  along  the  rod  and  at  the  boundaries.  These  first  two  examples 
also  show  clearly  the  nature  of  the  couphng  introduced  by  the  moderate  rotation 
theory  as  compared  to  the  infinitesimal  theory,  which  is  primarily  that  an  extensional 
response  is  induced  by  any  flexural  motions.  In  the  last  example  with  moderate 
rotation,  we  showed  that  the  theory  is  not  valid  for  vibration  in  the  fifth  flexural 
mode  even  though  it  was  valid  for  vibration  in  the  first  flexural  mode  for  the  same 
rod.  This  highlighted  the  fact  that  the  theory  is  obviously  not  valid  for  arbitrary 
deformations  of  the  rod. 

Our  last  example  was  the  whirhng  rod.  We  developed  the  balance  laws  for  a  rod 
rotating  about  an  end  at  a  constant  angular  velocity.  Our  intention  was  to  provide 
a  firm  foundation  for  the  problem,  which  has  a  long  but  inconsistent  history.  We 
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established  balance  laws  for  a  large,  steady  modified  auxiliary  motion  and  for  any 
superposed  vibrations.  We  did  not  solve  for  the  large  motion  or  for  the  vibration 
modes  because  we  left  the  explicit  form  of  the  free  energy  to  be  arbitrary.  However, 
we  discussed  an  approximate  solution  procedure  that  can  be  used  for  this  purpose. 
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